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FOREWORD 


This  report  documents  research  into  novel  uninhabited  aerial  system  reconfigurable  structures. 
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1  EXECUTIVE  SUMMARY 

Reconfigurable  structures  have  been  the  topic  of  much  research  in  recent  years.  Articulation, 
mechanisms,  skin  materials,  and  sensory  systems  to  increase  maneuverability  and  flight  mission 
performance  have  all  been  studied.  In  this  report,  a  particular  type  of  artificial  hair  sensor  is 
studied  for  use  as  a  flow  sensor  for  advanced  maneuvering.  Additionally,  skins  for 
reconfigurable  wings  are  considered  in  terms  of  a  cellular  structure  filled  with  a  variable  stiffness 
material,  such  as  honeycomb  filled  with  shape  memory  polymer  (SMP).  The  presented  work 
outlines  an  investigative  effort  to  determine  the  feasibility  of  such  a  system,  mechanically 
characterize  the  system  given  chosen  candidate  materials,  and  develop  an  optimization  scheme 
for  tailoring  the  system  to  specific  applications.  The  chosen  candidate  materials  of  epoxy  SMP 
and  aluminum  honeycomb  are  individually  characterized.  Using  these  materials,  the  Young’s 
moduli  and  shear  moduli  of  the  filled  honeycomb  composite  are  also  experimentally  determined 
above  and  below  the  transition  temperature  of  the  polymer.  These  results  are  then  used  to 
validate  an  analytical  model  taking  into  account  the  geometry  of  the  honeycomb.  A  finite 
element  analysis  (FEA)  is  also  presented  for  comparison  with  the  experimental  results  and 
analytic  model.  Once  validated,  the  analytic  model  is  then  used  in  an  optimization  scheme  to 
optimize  the  honeycomb  geometry  given  an  assumed  set  of  constraints  and  desired  skin 
properties.  Finally,  limitations  of  each  element  of  the  analysis  and  future  work  are  discussed. 
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2  ENERGY-BASED  RECONFIGURABLE  MICRO-AERIAL  VEHICLE  (MAV) 

FLIGHT  STRUCTURES  DESIGN 

2.1  Background 

Although  robust  control  strategies  have  improved  flight  system  performance  in  the  presence  of 
unmodeled  disturbances  such  as  gusts,  aircraft  are  nonetheless  susceptible  to  disturbances  that 
cause  significant  departure  from  linear  behavior.  In  this  work,  artificial  hair  sensor  array 
feedback  is  combined  and  utilized  in  a  recurrent  neural  network  (RNN)  approach  to  predict 
airfoil  loading  as  a  consequence  of  incident  gusts.  It  is  demonstrated  that  with  proper  recurrent 
neural  network  architecture,  distributed  sensory  array  information  can  be  effectively  combined 
and  utilized  to  predict  disturbances  in  lift  and  moment.  Such  results  indicate  that  empirical 
learning-based  approaches  to  utilizing  uncharacterized,  massive  arrays  of  flow  sensors  may  be 
both  practical  and  beneficial. 

2.2  Introduction 

Despite  their  success,  engineered  flight  systems  have  yet  to  achieve  the  performance  of 
biological  flying  systems  in  terms  of  maneuverability,  efficiency  and  robustness  to  disturbances. 
Robust  control  strategies  for  flight  systems  have  been  developed  that  provide  improved  stability 
and  performance;  however,  such  systems  are  susceptible  to  unmodeled  disturbances  such  as 
gusts.  In  this  work,  an  approach  to  disturbance  awareness  using  RNNs  and  distributed  artificial 
hair  sensor  feedback  is  developed.  Specifically,  an  RNN  is  used  to  estimate  body  reaction  (in 
terms  of  force  and  moment)  to  flow  disturbances  using  spatially  distributed  artificial  hair  sensor 
feedback. 

Biological  systems  present  a  variety  of  sensory  adaptations  that  potentially  provide  the  same,  or 
very  similar,  information  engineered  flight  systems  require  to  achieve  performance  similar  to 
their  biological  counterparts.  For  instance,  the  bat  may  use  muscle  and  hair  sensor  feedback  to 
achieve  extraordinary  flight  performance  [1,  2].  In  fact,  hair  sensors  may  be  the  most  sensitive 
sensory  modality  developed  in  nature  [3].  In  response  to  these  observations,  engineers  have 
begun  developing  insect-grade  flow  sensors  -  cheap,  numerous  and  uncharacterized  sensors  - 
and  investigating  potential  benefits  of  fly-by-feel  flight  systems.  Algorithms  and  hardware  that 
can  deal  effectively  with  massively  distributed  and  uncharacterized  sensory  arrays  will  be 
necessary  for  such  systems  to  come  to  fhiition. 

Engineered  flying  systems  currently  utilize  sensors  for  situational  awareness  and  state/parameter 
estimation,  but  the  number  of  sensors  is  limited  and  the  sensors  are  comparatively  expensive.  As 
such,  emphasis  is  placed  on  point-wise  sensor  accuracy  rather  than  redundancy.  Sensor 
redundancy  is  not  only  useful  but  necessary  for  biological  systems  [4,  5].  There  are  several 
reasons  to  consider  highly  redundant,  distributed  sensory  systems  for  flight  control.  First  and 
foremost,  redundancy  provides  robustness  to  sensor  failure  -  a  systemic  problem  in  systems  with 
very  few  high-grade  sensors.  Redundancy  also  provides  more  information  regarding  estimated 
quantities.  In  contrast  to  engineered  systems,  biology  generally  employs  enormous  levels  of 
redundancy  in  both  control  and  estimation  systems  [6]. 

As  described  later,  the  flow  features  of  interest  here  take  the  form  of  wind- vector  fields  over  an 
aircraft.  Such  features  are  inherently  nonlinear,  spatially  distributed,  and  unsteady.  An  accurate 
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physics-based  model  that  describes  these  features  exists,  but  a  real-time  solution  of  the  nonlinear 
Navier-Stokes  equations  is  computationally  intractable  at  this  time.  Nevertheless,  such  flow 
features  are  easily  observable  by  a  variety  of  sensor  types  and  lend  themselves  immediately  to 
very  successful  empirical  approaches  that  utilize  them  directly  for  flow  and  aircraft  state 
estimation  [7,  8].  Flow  fields  incident  on  an  aircraft  contain  information  regarding  aircraft  state 
and  relative  wind.  Temporal  variations  in  flow  topology  may  also  serve  as  a  means  for  prediction 
of  fluid-body  interaction.  For  instance,  temporal  and  spatial  variations  in  distributed  hair  sensor 
signals  may  provide  reliable  prediction  of  an  impending  consequential  flow-state  change.  At  the 
very  least,  such  signals  are  intimately  related  to  aircraft  body  forces  and  moments,  as  described 
in  [7,  8]  and  the  references  therein.  Hence,  RNNs  are  investigated  as  a  means  to  provide  both 
state  and  parameter  estimates,  as  well  as  a  means  of  predicting  future  body  forces  and  moments  - 
a  capability  conducive  to  a  fly-by-feel  paradigm. 

2.3  Flow  Features,  Artificial  Hair  Sensors,  and  Distributed  Sensor  Arrays 

All  flow  sensing  hairs  have  been  found  on  many  biological  models,  including,  but  not  limited  to, 
crickets,  spiders,  locusts,  fish,  and  bats  [1,9,  10].  Such  sensory  adaptations  are  used  for  predator 
evasion,  environmental  awareness,  and  flight  performance.  Even  blind  fish  have  demonstrated 
the  ability  to  school  by  using  distributed  hair-like  sensors,  a  behavioral  fact  that  indicates  the  use 
of  distributed  flow  measurements  [11].  Like  most  biological  sensory  modalities,  these  systems 
are  generally  fairly  robust  due  to  redundancy,  material  properties,  directional  sensitivity,  and  the 
adaptive  nature  of  neurological  systems.  Hair  sensors  are  a  fairly  new  addition  to  the  suite  of 
anemometers,  and  function  by  measuring  the  mechanical  response  of  a  beam-like  extension  into 
a  flow  field.  Drag  from  hair-flow  interaction  results  in  forces  and  moments,  measured  at  the  base 
of  the  hair,  which  are  a  function  of  the  local  flow  velocity  field.  Output  from  hair  sensors  is  a 
result  of  drag  along  the  fully  wetted  hair  surface  and  can  be  calculated  as  the  moment  integral  of 
drag  on  the  hair.  A  noteworthy  review  of  hair-like  sensors  is  provided  in  [10],  with  interesting 
experimental  and  theoretical  analysis  provided  by  [12,  13]  and  the  references  therein.  Some 
groups  have  devised  several  designs  for  single  and  arrays  of  artificial  hair  flow-sensors,  some  of 
which  have  been  implemented  successfully  in  flow  experiments  [10,  14].  Until  recently,  most  of 
the  research  has  been  limited  to  material  and  mechanical  design  studies.  Several  groups  have 
now  begun  investigating  identification  of  aerodynamic  phenomena  and  aerodynamic 
control/estimation  strategies  that  utilize  artificial  hair  sensors  [15,  16].  Computational  analysis 
has  yielded  promising  results  indicating  that  such  sensors  can  be  used  for  boundary  layer  and 
flow  estimation.  In  this  work,  we  provide  preliminary  results  indicating  AHS  do  offer  the 
capability  for  flow  estimation  and  flow  feature  identification. 

Furthermore,  hair  sensors  offer  the  potential  for  observation  beyond  point  measurements, 
providing  direct  access  to  observable  flow  features  that  result  from  fluid-structure  interaction. 
Intuition  and  research  suggest  that  flying  animals  likely  react  to  flow  structures  as  they  are 
encountered  and  are  doing  so  based  on  previous  experience  and  memory.  Such  flow  structures 
may  be  identifiable  by  observing  patterns  and  features  in  the  flow  velocity  field  near  and  upon  a 
body  in  flight.  Thus,  hair  sensors  may  directly  sense  flow  structures  (and  temporal  variations 
thereof)  as  important  features,  or  observables,  to  provide  an  organism  a  predictive  capability  to 
adjust  control  strategies  to  achieve  a  goal  even  in  the  face  of  changing  flow  conditions  and 
disturbances. 


3 

Approved  for  public  release;  distribution  unlimited. 


When  distributed  arrays  of  flow  sensors  are  used,  observation  of  such  features  becomes  similar 
to  image  or  pattern  recognition  problems  of  artificial  intelligence.  Spatially  distributed  sensors 
enable  imaging  of  flow  feature  observables,  including  stagnation,  separation,  and  reattachment 
points  (lines)  such  as  those  illustrated  in  Figure  1 .  Such  features  are  intimately  related  to  flight 
body  state  and  stability  and  can  thus  be  used  for  both  control  and  estimation  of  meaningful 
parameters  and  states  [8],  For  instance,  Mangalam  et  al.  have  developed  a  strategy  for  direct 
observation  and  real-time  extraction  of  critical  aerodynamic  flow  characteristics,  loads,  and 
moments  from  surface  signatures  [7,  8],  In  their  approach,  hot-film  anemometry  is  used  to 
directly  measure  location  of  the  flow  features  mentioned  above,  along  with  flow  transition, 
turbulence,  shock,  vortices,  etc.  Related  approaches  such  as  that  by  Goman  and  Khrabrov  in  [17] 
demonstrated  that  unsteady  lift,  drag,  and  moments  can  be  expressed  as  functions  of  the 
instantaneous  separation  point  location,  angle  of  attack,  and  pitch  rate. 


Figure  1:  Stagnation,  Separation,  and  Reattachment  Flow  Features 

Fundamentally,  all  practical  approaches  can  be  described  as  striving  to  introduce  a  priori 
knowledge  into  an  estimation  scheme  that  will  uniquely  and  robustly  determine  useful 
aerodynamic  parameters.  Therefore,  it  seems  logical  to  consider  whether  or  not  the  problem  of 
aircraft  state  estimation  and  control  can  be  enhanced  by  using  a  nonlinear  system  identification 
approach  that  utilizes  a  plethora  of  redundant  sensors,  rather  than  relying  upon  complex 
mathematical  models  requiring  online  numerical  integration.  After  all,  it  seems  nature  is  already 
solving  the  nonlinear  equations  of  motion  exactly;  why  burden  engineered  flight  systems  with 
mere  approximations  when  such  systems  can  instead  feel  their  way  through  flight  and  at  times, 
learn  and  adapt  through  experience.  The  complexity  of  dealing  with  hundreds,  if  not  thousands, 
of  spatially  distributed  sensors  is  traditionally  either  computationally  burdensome  or  relies  upon 
extensive  offline  statistical  analysis  and  system  identification.  Although  adaptive  estimation  and 
control  methodologies  exist  for  such  applications,  these  generally  still  require  physics-based 
models  that  do  not  exist  (and  may  never)  for  such  nonlinear  problems.  Neural  networks  stand  as 
a  strong  candidate  for  alleviating  the  above  challenges. 

Consider  a  large  array  of  sensors,  with  distributed  neural  networks  to  combine  and  adapt  for  each 
local  group  of  sensors,  the  signals  of  which  feed  upward  through  an  architecture  similar  to  that 
shown  in  Figure  2.  With  sufficient  theoretical  understanding,  the  majority  of  computation  could 
be  performed  locally.  Local  computation  results  could  be  fed  upward  into  a  traditional  control 
design,  modified  to  utilize  such  data.  It  is  the  author’s  opinion  that  traditional  control 
methodologies  shouldn’t  be  replaced,  but  rather  augmented  with  rich  sensory  feedback  and 
learning-based  architectures  such  as  that  presented  here. 
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Figure  2:  Biologically  Inspired  Sensory  Architecture 
2.4  Artificial  Hair  Sensor  Principles 

Artificial  hair  is  strongly  biologically  inspired.  It  is  now  well-known  that  biological  models 
utilize  wind  chill  and/or  microscopic  hairs  for  wind  orientation,  turbulence  detection,  among 
other  adaptations.  For  example  bats  may  use  hair  for  flow  detection  applications.  D’Angelo  et  al. 
investigated  bats’  use  of  hair  for  turbulence  and  consequences  of  hair  removal  upon  bat  flight. 
Bat  flight  metrics  experience  degradation  upon  hair  removal  [1].  Other  animals  utilizing  hair  or 
hair-like  structures  for  sensing  flows  include  insects  and  fish  [12,  11].  Hair  as  a  sensor  is  a 
successful  strategy  in  biology  for  multiple  reasons  including  low  mass-density  and  regrowth 
capability. 

In  light  of  the  aforementioned  biological  discoveries,  our  group  has  designed  and  constructed 
artificial  hair  sensors  from  microscale  glass  fibers  and  carbon  nanotubes.  Figure  3  depicts  a  hair 
sensor  undergoing  deflection  due  to  drag  forces.  As  the  hair  is  effectively  cantilever  from  a 
surface,  it  protrudes  within  and  potentially  out  of  a  boundary  layer,  as  illustrated  alongside  the 
hair.  Thus,  fluid  motion  induces  drag  upon  the  hair,  causing  a  moment  and  shear  force  to  exist  at 
the  hair  base.  By  growing  carbon  nanotubes  along  the  hair  contained  within  the  glass  pore,  hair 
deflections  can  be  made  to  cause  changes  in  resistance  due  to  nanotube  deformations. 
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Flow  Velocity 
Distribution 


Figure  3:  Artificial  Hair  Function  Concept 

To  be  specific,  if  a  coordinate  axis  is  set  up  such  that  elevation  along  the  hair  is  represented  by  z, 
total  hair  height  as  and  drag  upon  the  hair  (drag  normal  to  the  hair  axis  as  a  function  of  z) 
represented  by  5(z),  the  total  drag  upon  the  hair  can  be  integrated  to  yield  (for  one  dimension) 

My  =  J  S(z)zdz  (1) 

0 

Zh 

Fx=  j  S(z)  dz  (2) 

0 

where  My,  are  the  moment  about  the  y  axis  and  force  along  the  x-axis,  respectively.  The  drag 
5  is  a  function  of  Reynolds  number  and  is  classically  defined  as 

5  =  ^pv^CdA  (3) 

where  p  is  fluid  mass  density,  v  is  the  magnitude  of  the  flow  velocity  vector  normal  to  the  hair 
center  axis  (scalar),  is  the  coefficient  of  drag  and  is  empirically  determined  to  combine  form 
and  skin  drag,  and  A  is  reference  area.  Thus,  by  measuring  moment,  force,  or  some  combination 
of  the  two,  the  velocity  of  flow  at  hair  location  can  be  estimated. 

2.5  RNNs  and  Real  Time  Recurrent  Learning 

Repeated  (in  space  or  time)  measurements  of  a  quantity  generally  provides  more  information 
regarding  the  quantity  of  interest.  In  general,  Bayes’  theorem  provides  a  systematic  approach  to 
solving  most  estimation  problems.  In  the  case  of  linear  systems  with  characterized  Gaussian 
noise  in  process  and  measurement,  the  well-known  Kalman  filter  provides  an  optimal  solution. 
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extracting  precisely  the  available  information  from  all  measurements.  However,  such  optimal 
closed-form  solutions  only  exist  in  very  constrained  circumstances.  Furthermore,  the  Kalman 
filter  is  also  problematic  in  high  dimensional  settings,  exhibiting  both  numerical  instability  and 
undesirable  computational  burden.  Nonlinear  problems  such  as  those  involving  flight  and  flow 
dynamics  generally  necessitate  development  of  approximate  estimators  or  learning-based 
approaches,  especially  those  for  which  accurate  models  do  not  exist  [18]. 

Making  a  difficult  problem  even  more  so,  a  disturbance  in  a  flow  can  occur  over  a  short  period 
of  time.  The  resulting  dynamic  prediction  problem  presents  a  difficulty  for  classical  disturbance 
observation  or  prediction  frameworks  (based  on  dual  or  joint  Kalman  filters,  etc.).  Such 
approaches  require  great  computational  speed  to  accurately  estimate  the  undesired  input  to  the 
flight  system  all  the  while  providing  sufficient  time  for  an  actuator  system  to  reject  it.  Systems 
such  as  light  detection  and  ranging  (LIDAR)  attempt  to  overcome  such  issues  by  directly 
observing  a  disturbance  ahead  of  the  craft  [19].  Control  strategies  can  then  be  used  to  mitigate 
the  impending  disturbance.  However,  LIDAR  systems  are  complex  and  costly.  The  ability  to 
predict  impending  disturbances  based  on  changing  local  flow  features  could  be  beneficial  to 
smaller  flight  systems  that  cannot  support  systems  like  LIDAR. 

The  success  of  biological  models  is  likely  due  to  the  utilization  of  prior  knowledge  in  the  form  of 
predictions  based  on  previous  experience  and  training.  Prior  knowledge  is  important  for  accurate 
estimation  and  prediction  in  complex  dynamic  environments  and  can  be  obtained  through 
machine  learning  techniques.  Furthermore,  in  contrast  to  model  based  control  and  estimation 
strategies,  the  prediction  of  external  disturbances  could  provide  a  solution  much  like  that  of 
LIDAR,  except  as  an  inferred  prediction. 

Considering  the  potential  complexity  of  flow  features,  massively  distributed  flow  sensor  arrays 
and  the  naturally  ill-posed  problem  of  dynamic  prediction,  machine-learning  techniques  and 
specifically  RNNs  appear  as  strong  candidates  for  a  practical  solution.  RNNs  are  fully  connected 
networks  containing  internal  feedback  and  time  delay,  extending  the  concept  of  standard  neural 
networks  to  nonlinear  dynamic  systems.  Figure  4  illustrates  a  RNN  consisting  of  an  input  layer, 
hidden  inter-connected  neurons,  and  output  layer.  However,  unlike  standard  neural  networks,  the 
hidden  units  are  fully  connected.  RNNs  can  provide  universal  approximation  for  any  nonlinear 
function  of  both  space  and  time,  provided  sufficient  hidden  neurons  [18]. 


Input  Layer  Hidden  Neurons  Output  Layer 


Figure  4:  Example  RNN  Architecture 
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Consider  a  noise-free  input  and  RNN  system  like  the  one  previously  illustrated.  Consider  a 
connection  between  any  two  neurons  at  time  step  n  with  states  and  respectively,  where  x^ 
is  a  vector  of  size  q  x  1.  The  weight  describing  the  connection  for  a  signal  propagating  from 
neuron  j  to  neuron  /  can  be  written  .  Consider  any  subset  of  neurons  in  a  network  as  being 
input  and  output  units,  with  state  vectors  denoted  as  and  y„  and  sizes  m  x  1  and  p  x  1, 
respectively.  Then,  the  operation  of  any  RNN  can  be  written  in  the  form  (directly  from  [18]). 

Xn+l  =  ^i^aXn  +  (4) 

Tn  =  ^cX-n  (5) 

where  Vk^is  a  p  X  q  input  weight  matrix,  and  0:R'^  ^  R'^is  the  activation  function  for  hidden 
units.  The  network  output  is  available  by  sampling  whichever  neurons  are  chosen  as  output 
nodes.  The  matrices  Vk^and  M4determine  the  connectivity  of  the  neurons,  which  can  include 
connectivity  between  inputs,  outputs,  hidden  units,  or  any  combination  thereof  The  matrix 
specifies  the  output  neurons  and  can  be  considered  a  measurement  model  in  the  sense  of  MIMO 
systems.  It  is  generally  populated  sparsely  by  ones,  which  return  output  layer  neuron  values. 

Here,  <I)  is  chosen  as  the  hyperbolic  tangent  for  hidden  units  and  the  linear  output  function  for  all 
output  layer  neurons.  Delays  are  implicit  in  the  structure  of  Equation  (4),  as  input  signals  to  the 
system  may  take  some  time  to  elicit  a  response  in  distant  neurons,  depending  on  the  sparsity  of 
the  matrix  W.  Also,  the  generality  of  the  above  form  allows  for  customization  of  network 
architecture,  trimming  of  specific  connections  (setting  relevant  entries  of  W  to  zero),  etc. 
Equation  (4)  takes  a  general  form  similar  to  discretized  equations  of  motion,  which  not  only 
intuitively  supports  the  aforementioned  universal  approximation  theorem,  but  suggests  why  it  is 
important  to  use  a  sufficient  number  of  hidden  neurons  to  approximate  a  given  system. 

2.6  Neural  Network  Training  Using  the  Extended  Kalman  Filter 

Many  algorithms  exist  for  training  RNNs.  Recently,  supervised  training  using  nonlinear 
sequential  state  estimators  has  emerged  as  one  of  the  most  successful  approaches  [18,  20].  Here, 
an  extended  Kalman  filter  (EKE)  is  used  to  train  the  connection  weights  in  a 
prediction/correction  strategy,  using  a  training  set  consisting  of  hair  sensor  output  and  desired 
network  response  (lift  and  moment  on  an  airfoil).  Because  inputs  and  desired  outputs  are  known 
a  priori,  the  network  becomes  a  function  of  the  weights  that  are  fully  observable  to  the  algorithm. 
Thus,  application  of  the  EKE  to  estimation  of  the  weights,  taken  here  as  system  states,  is 
straightforward,  but  does  require  computation  of  the  gradient  of  output  with  respect  to  weights 
[18].  From  Haykin  [18],  the  algorithm  can  be  summarized  as  follows.  Given  the  training  sample 

set  f  « ’  « >n-\  ^  where  u„  is  the  input  vector  applied  to  the  RNN  and  d„  is  the  corresponding 
desired  response,  direct  application  of  the  EKF  to  the  dynamic  system  composing  the  RNN 
becomes: 


+  Qy,n  ^  (6) 

(7) 
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where  G„  denotes  the  Kalman  gain  matrix,  B  denotes  a  linearized  measurement  matrix,  b„  is  a 
function  denoting  the  RNN,  w  is  a  vector  composed  of  RNN  connection  weights,  P  is  the 
measurement  error  covariance  matrix,  and  is  a  measurement/process  noise  covariance 
matrix.  The  reader  interested  in  a  detailed  derivation  and  discussion  of  the  above  algorithm 
should  consult  Haykin  [18].  In  the  next  section,  a  potential  flow  gust  simulation  is  briefly 
described  that  provides  approximate  hair  sensor  output  alongside  airfoil  force  and  moment  time 
history.  Taken  together,  this  information  forms  a  set  of  training  data  for  a  RNN. 


2.7  Quasi-Steady  Potential  Flow  Simulation 

In  order  to  approximate  a  flow  disturbance,  such  as  a  gust,  airfoil  loading  due  to  nearby 
convection  of  vortex  structures  in  a  potential  flow  is  simulated.  The  approach  used  here  is  largely 
based  on  theory  in  [21],  is  similar  to  the  approach  found  in  [22],  and  consists  of  a  Joukowski 
airfoil  in  a  flow.  A  pair  of  vortices  is  introduced  into  the  flow,  one  contained  within  the  airfoil 
and  one  external.  The  position  of  the  internal  vortex  is  calculated  as  to  ensure  the  airfoil 
boundary  is  a  streamline,  effectively  satisfying  a  non-penetrative  wall  boundary  condition.  The 
vortex  propagates  downstream  according  to  Routh’s  rule  [23],  passing  by  the  airfoil.  Total  airfoil 
circulation  is  modified  at  each  time  step  such  that  the  Kutta  condition  is  satisfied,  resulting  in 
loading  of  the  airfoil.  A  second  order  Runge-Kutta  integration  scheme  is  used  to  compute  a 
solution  to  the  vortex  convection.  During  execution  of  this  algorithm,  local  flow  velocity  along 
the  airfoil  is  computed  at  discrete  locations,  chosen  a  priori  as  artificial  hair  sensor  locations.  The 
magnitude  of  these  velocities  is  then  computed,  emulating  response  of  an  array  of  directionally 
insensitive  artificial  hair  sensors.  Consider  a  cylinder  of  radius  a,  located  at  point  ^  in  an 
irrotational,  two-dimensional  flow  of  incompressible  inviscid  fluid  in  what  will  be  referred  to  as 
the  ^-plane.  Then,  by  potential  flow  theory,  the  complex  potential  describing  the  i/-plane  flow 
can  be  written  as 


w  =  f/ooe‘“((  -  (o)  +  f/ooa2e-‘«/((  -  ^o)  (H) 

where  a  is  the  angle  of  incidence  of  the  free  stream  with  respect  to  the  horizontal  axis.  If  a 
vortex  of  strength  k  is  located  at  and  the  cylinder  has  circulation  of  strength  Kq,  then  the 
potential  becomes 


UooCi^e 

w  =  f/ooe'“«'  -  ^o)  +  +  iKo  log«'  -  ^0  -  ^i)  -  i^ilogC^  -  Co  “  ■^)  (12) 

with  the  additional  terms  representing  the  contribution  of  the  cylinder  circulation,  the  vortex  at 
point  ^1,  and  a  matching  vortex  inside  the  cylinder  ensuring  the  cylinder  boundary  is  a 
streamline.  If  the  external  vortex  is  allowed  to  convect  according  to  the  external  flow,  the  terms 
in  Equation  (12)  must  be  modified  according  to  a  numerical  integration  routine.  To  compute  the 
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flow  about  a  Joukowski  airfoil,  a  coordinate  transformation  is  applied,  resulting  in  a  conformal 
mapping.  This  mapping  is  generally  called  Joukowski ’s  transformation,  and  takes  the  form 


z  = 


Computation  of  the  velocity  field  about  the  airfoil  (in  the  z-plane)  is  straightforward,  and 
immediately  results  from  differentiation  of  Equation  (12)  with  respect  to  z.  From  the  resulting 
transformed  flow,  the  velocity,  pressure,  and  therefore  forces  and  moments  acting  upon  the 
airfoil  can  be  computed.  Hair  sensor  response  is  taken  to  be  completely  in  phase  and  proportional 
to  the  magnitude  of  flow  velocity  at  discrete  points  along  the  airfoil  boundary.  The  theorem  of 
Blasius  [21]  is  used  to  calculate  forces  and  moments  on  the  airfoil.  Figure  5  depicts  a  snapshot  of 
a  convecting  vortex  near  an  airfoil.  The  imaginary  component  of  the  flow  potential  is  plotted  as  a 
contour  surface,  the  contours  of  which  describe  streamlines.  Note  that  although  the  vortex  has 
severely  modified  the  flow,  the  airfoil  boundary  is  a  streamline.An  example  result  of  vortex 
convection  near  an  airfoil  is  illustrated  in  Figure  5  and  Figure  6,  in  which  an  airfoil  at  zero  angle 
of  attack,  Reynolds  number  of  3.5  x  10^,  is  exposed  to  a  vortex  as  it  convects  downstream.  The 
initial  position  of  the  vortex  was  chosen  by  sampling  a  uniform  distribution,  modified  to  ensure 
that  the  vortex  would  not  convect  closer  than  five  chord  lengths  to  the  airfoil.  The  initial  distance 
upstream  is  specified  as  100  chord  lengths.  Strength  of  the  vortex  was  also  generated  randomly 
to  have  vorticity  -10  s'^  <  ki  <  10  s'\ 


-6  -4  -2  0  2  4  6 


Y 

Figure  5:  Snapshot  of  Vortex  with  Airfoil  (Imaginary  Component  of  Potential) 

The  lift  and  moment  resulting  from  such  a  vortex  gust  is  illustrated  in  Figure  5,  with  leading 
edge  hair  sensor  responses  illustrated  in  Figure  6  and  Figure  7  as  an  interpolated  surface.  One 
can  readily  identify  movement  of  forward  stagnation  point  in  Figure  7  as  a  dark  region.  The 
movement  of  stagnation  point  is  related  to  airfoil  loading,  as  can  be  inferred  visually  by 
comparing  the  major  features  of  Figure  5,  Figure  6,  and  Figure  7.  Only  a  small  time  span  has 
been  plotted  which  contains  the  relevant  information.  Fift,  moment  and  hair  sensor  output  all 
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asymptotically  approach  zero  as  one  moves  away  from  the  time  of  maximum  airfoil-vortex 
interaction. 


Figure  6:  Lift  and  Moment  Time  History 
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Figure  7:  Hair  Sensor  Response  vs.  Time 


2.8  Simulation  Results 

The  potential  flow  simulation  described  previously  was  then  used  to  produce  training  sets  by 
randomly  selecting  vortex  strength  and  initial  position,  but  with  airfoil  angle  of  attack  maintained 
at  zero  (to  simplify  preliminary  study)  and  free-stream  velocity  dictated  by  a  specified  Reynolds 
number.  As  such,  parameters  were  non-dimensionalized  according  to  chord  length,  and  the 
Reynolds  number  was  set  at  3.5106.  For  the  present  investigation,  hair  length  is  assumed  to  be 
sufficient  to  be  dominated  by  flow  external  to  any  boundary  layer  that  would  actually  be  present. 
Obviously,  a  hair  fully  contained  within  a  boundary  layer  can  no  longer  be  assumed  to  provide 
such  a  response,  but  the  present  learning-based  approach  is  largely  an  effort  to  render  such 
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modeling  complexity  irrelevant.  A  training  set  serves  to  manifest  a  physically  meaningful 
relationship  between  hair  sensor  12  response  and  aerodynamic  loading.  The  goal  here  was  to 
construct  a  nonlinear  model  (a  RNN)  that  approximates  the  complex  hair  sensor-flow 
relationship  with  regard  to  external  disturbances,  while  encouraging  a  predictive  capability.  A 
RNN  architecture  was  chosen  which  consists  of  an  input  layer  of  size  (n  +  3),  where  n  is  the 
number  of  hair  sensors,  with  three  additional  input  neurons  for  bias  and  net  output  feedback  (lift 
and  moment  estimates).  For  the  results  presented  here,  three  hair  sensors  were  used,  extending 
from  mid-chord  on  both  upper  and  lower  airfoil  surfaces,  equally  spaced  around  the  leading 
edge.  The  RNN,  consisting  of  six  input  neurons,  six  hidden  neurons  and  two  output  neurons,  was 
trained  on  the  aforementioned  gust  simulation  data.  Preliminary  results  indicate  success  as  far  as 
training  a  RNN  to  accurately  describe  airfoil  loading  as  a  function  of  hair  sensor  feedback. 

Figure  8,  Figure  9,  and  Figure  10  depict  training  results  from  a  single  gust  simulation.  Using  the 
hair  output  of  Figure  8,  a  RNN  can  be  trained  by  repeatedly  looping  through  the  time  history  of 
the  training  set.  Here,  repeated  loops  are  referred  to  as  epochs.  The  total  error  in  time  (Euclidean 
norm)  is  illustrated  for  each  epoch  in  Figure  9.  Within  80  epochs,  the  system  converged  to  a 
degree  resulting  in  imperceptible  difference  between  RNN  and  desired  response.  While  such  a 
training  methodology  would  likely  yield  poor  generalization  to  other  gusts,  it  demonstrates  the 
effectiveness  of  the  training  algorithm  chosen. 


Figure  8:  Hair  Output  Time  History 
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Figure  9:  Training  Error  vs.  Epoch 

The  desired  output  of  the  RNN  as  a  function  of  hair  sensor  input  agrees  well  with  the  desired  lift 
and  moment  functions,  as  can  be  seen  in  Figure  10.  The  real  question  is  whether  or  not  such  a 
network  can  predict  the  lift  and  moment  for  a  general  vortex.  Such  questions  will  be  addressed  in 
the  final  manuscript. 


Figure  10:  Desired  Output  vs.  Network  Output  After  Training 

We  anticipate  results  that  will  describe  a  network  architecture  that  uses  distributed  hair  sensor 
signals  to  predict  the  forces  and  moments  resulting  from  a  gust.  Such  an  architecture  will  utilize 
a  two-component  neural  network  system,  the  first  consisting  of  hair  sensor  fusion  layers  which 
are  not  recurrent,  but  provide  an  internal  output  to  a  RNN  that  is  dedicated  to  prediction.  Such  an 
architecture  serves  as  a  foundation  for  investigating  distributed  sensor  compression  and 
associative  memory  design  for  discovering  critical  flow  features,  such  as  stagnation  point,  and 
for  modeling  in  an  effective  and  practical  manner  the  complex  physics  of  the  problem.  The 
proposed  architecture  may  also  lend  itself  to  incorporating  other  sensory  modalities  such  as 
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pressure  and/or  structural.  As  the  potential  flow  model  used  in  this  work  was  fairly  crude,  a  high- 
fidelity  Navier-Stokes  simulation  will  be  used  to  generate  equivalent,  but  much  more  realistic 
hair  sensor  and  airfoil  loading  data  for  future  investigation.  Such  a  simulation  includes  unsteady, 
complex  flow  phenomena  such  as  turbulence  and  separation.  Indeed,  it  is  these  phenomena  that 
are  of  the  utmost  interest  for  this  work,  and  the  results  that  come  from  such  a  study  will  be 
extremely  interesting. 

2.9  Experimental  System  Demonstration 

To  demonstrate  the  functionality  of  the  hair  sensors  and  the  application  of  the  above  neural 
network  control  algorithms,  the  flow  characteristics  around  a  cylinder  in  cross-flow  was 
experimentally  measured  using  a  hair  sensor.  A  single  hair  sensor  was  mounted  to  a  cylinder 
and  the  resistance  across  the  top  and  bottom  electrode  measured  at  various  angles  as  the  cylinder 
was  rotated  in  the  flow.  Sensor  output  at  each  angle  measured  is  shown  in  Figure  11.  Relatively 
stable  signal  was  achieved  at  each  angle  as  shown  by  the  limited  amount  of  variation  along  the 
Time  axis  and  the  standard  deviation  for  each  angle  shown  in  Figure  12. 


Figure  11:  Sensor  Response  Surface 
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Figure  12:  Sensor  Standard  Deviation 

Relative  flow  velocity  around  the  cylinder  using  PIV  is  shown  in  Figure  13.  Interior  to  the 
cylinder  is  the  output  from  the  hair  sensor.  Flow  velocity,  the  stagnation  point,  and  separation 
point  are  accurately  measured. 


X 


Figure  13:  Experimental  PIV  Wind  Tunnel  Data 
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2.10  Artificial  Hair  Sensor  Conclusion 

While  preliminary,  the  main  conclusion  that  can  be  garnered  from  this  work  is  that  RNNs  and 
EKF -based  training  are  very  effective  for  learning  overall  system  behavior.  A  challenge, 
however,  lies  in  generating  a  successful  network  architecture  to  achieve  a  desired  end.  Indeed, 
this  is  one  of  the  main  problems  with  using  neural  networks  in  practical  application.  However,  as 
was  demonstrated,  even  a  simple  approach  can  nonetheless  yield  valuable  results  for  a  problem 
that  might  otherwise  be  intractable.  The  nonlinear  modeling  capabilities  of  RNNs  are 
extraordinary  and  have  begun  to  evoke  excitement  in  the  robotics  and  machine  learning 
community  due  to  their  learning  and  memory  potential.  However,  their  application  to  flight 
estimation  and  control  applications  is  largely  uninvestigated.  With  the  potential  to  offer  both 
predictive  and  accurate  nonlinear  system  approximation  capabilities,  RNNs  may  yield 
performance  and  robustness  enhancements  that  were  previously  unobtainable.  Future  work  will 
address  the  concepts  and  ideas  mentioned  throughout  this  document  with  the  end  goal  of  testing 
such  strategies  in  wind  tunnel  experimentation.  How  one  might  utilize  the  present  RNN  approach 
in  current  aircraft  control  systems  will  also  be  addressed.  The  success  of  the  current  control 
paradigm  is  unquestionable,  but  methods  of  incorporating  flow  features  into  such  a  framework 
have  been,  to  date,  entirely  ad  hoc.  As  such,  RNNs  may  offer  a  consistent  empirical  approach  to 
dealing  with  massive  amounts  of  complex,  dynamic  sensory  data  to  improve  aircraft 
performance. 
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3  RECONFIGURABLE  SKIN  SYSTEM  DESIGN 


3.1  Introduction 

Reconfigurable  structures  have  been  the  topic  of  much  research  in  recent  years  including 
candidate  materials  and  systems  used  to  cover  them,  aka  skins.  It  has  also  been  proposed  that 
one  such  skin  could  be  made  of  a  cellular  structure  fdled  with  a  variable  stiffness  material,  such 
as  honeycomb  fdled  with  shape  memory  polymer  (SMP).  [24,  25,  26,  27] 

3.2  Concept  Description  and  Motivation 

The  aerodynamic  performance  of  aircraft  for  specific  mission  profile  segments  can  be  improved 
through  changes  in  wing  shape  [28].  Andersen  showed  that  in  a  variety  of  mission  segments 
such  as  dash,  cruise,  climb,  and  loiter  a  different  wing  configuration  optimized  aerodynamic 
performance.  The  current  goal  is  to  create  a  monolithic  skin  that  does  not  wrinkle  when 
deformed  and  is  able  to  strain  more  than  traditional  materials  while  supporting  aerodynamic 
loads. 

Cellular  structures  can  be  designed  with  relatively  low  in-plane  and  high  out-of-plane  stiffness. 
Low  stiffness  is  obtained  by  the  empty  space  within  each  cell,  providing  space  for  the  nearby 
thin  structures  to  deform.  Deformation  of  these  thin  beam  structures  can  be  inhibited  by  filling 
the  empty  space  or  by  adhering  face  sheets.  Thus  the  stiffness  of  the  cellular  composite  can  be 
controlled  by  judicious  choice  of  materials.  If  a  phase  change  or  variable  stiffness  material  such 
as  shape  memory  polymer  is  used,  the  effective  stiffness  of  the  skin  is  controllable  both  in 
magnitude  and  direction.  By  manipulating  the  stiffness  of  individual  cells,  for  example  by 
thermal  means  such  as  shown  in  Figure  14,  in  specific  patterns  the  global  material  properties  of 
the  skin  can  be  tailored  to  the  specific  in-flight  needs  of  the  aircraft.  A  few  potential  patterns  are 
shown  in  Figure  15.  One  of  the  goals  of  the  presented  research  is  an  experimental  investigation 
into  the  possible  range  of  effective  skin  properties,  representing  bounded  limits  of  the  system. 

As  a  representative  system,  epoxy  shape  memory  polymer  and  aluminum  honeycomb  were  used 
to  demonstrate  the  system  and  validate  the  developed  models. 
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Figure  14:  Filled  Honeycomb  Concept 


Figure  15:  Conceptual  Heating  Patterns 

Thus  the  properties  of  the  proposed  filled  honeycomb  reconfigurable  skin  require  experimental 
validation,  modelling,  and  optimization.  This  was  accomplished  by  experimentally 
characterizing  the  epoxy  SMP,  aluminum  honeycomb,  and  SMP/honeycomb  composite  at 
ambient  and  above  Tg  temperatures  in  both  in-plane  cardinal  directions.  This  data  along  with 
FEA  was  then  used  to  validate  an  analytical  model.  Finally,  using  the  experimentally  validated 
analytical  model,  topology  optimization  on  the  composite  skin  was  performed  to  find  the  ideal 
honeycomb  geometry  given  a  prescribed  set  of  boundary  conditions. 

One  of  the  major  remaining  barriers  to  the  widespread  adoption  of  thermally  activated  shape 
memory  polymer  is  the  method  used  to  heat  them.  Presented  is  an  investigation  into  using  5  pm 
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nickel  particles  aligned  into  chains  as  embedded  Joule  heaters  for  epoxy  SMP.  The  high  density 
of  particle  chain  heaters  reduces  the  time  and  energy  required  to  reach  transition  by  minimizing 
excess  heat  required  due  to  the  low  thermal  conductivity  of  the  polymer  by  heating  the  material 
more  uniformly.  The  chains  are  formed  by  curing  the  polymer  in  a  uniform  magnetic  field 
generated  by  two  sets  of  N42SH  Neodymium  magnets  above  and  below  the  sample.  Modeling 
of  the  induced  magnetic  field  within  and  between  particles  during  curing  and  an  analytical  model 
predicting  particle  mobility  in  a  fluid  with  respect  to  vibration  frequency  and  amplitude  are 
presented  and  discussed  in  context  to  the  current  work.  Since  epoxy  resin  has  a  high  viscosity, 
particle  mobility  is  encouraged  by  sonicating  the  sample  at  300  Hz  at  an  amplitude  of 
approximately  50  pm  prior  to  polymerization  using  an  industrial  shaker  and  Teflon  guides. 
Copper  mesh  electrodes  are  attached  to  the  resulting  samples  using  10  percent  by  volume  nickel 
particle  SMP  epoxy.  Significant  particle  alignment  is  confirmed  via  optical  microscope  images. 
Electrical  resistivity  is  measured  as  low  as  57  Q-mm  at  Nickel  volume  concentrations  of  1.0 
percent.  Infrared  images  of  the  samples  during  heating  are  presented  and  electrical  energy 
required  with  respect  to  sample  thermal  capacity  estimated. 

The  Morphing  Aircraft  Structures  (MAS)  phase  II  program  is  used  to  demonstrate  the 
capabilities  of  the  system  and  as  a  realistic  design  optimization  criteria.  Program  requirements 
are  listed  in  Table  1  and  include  aerodynamic  load  specifications,  desired  skin  dimensions, 
weight  limits,  and  skin  functionality  requirements. 

Table  1:  MAS  Program  Skin  Requirements 


MAS  Requirements _  English 

Total  Skin  Weight  <  0.95  Ib/sq-fl 

Out-of-plane  Air  Loads  400  Ib/sq-ft 

Nominal  Panel  Size  15”  by  20” 

Max  Out-of-plane  Deflection  0.1” 

Shear  Capability  30-75  degrees 

No  wrinkling  of  skin  in  any  position 
Enforced  displacement  in-plane 
Aerodynamic  loads  out-of-plane 


SI 


<4638.31  g/m2 
19,152.1  Pa 
0.381  by  0.508  m 
2.54  mm 

~  p/4  -  2p/3  radians 


3.3  Embedded  Heating  Elements 

Reconfigurable  skins  have  been  proposed  for  a  multitude  of  applications  including  morphing 
aircraft,  deployable  satellites,  and  shape  changing  cars.  One  of  the  main  correlations  between 
these  applications  is  the  need  for  a  skin  material  with  controllable,  adaptable  material  properties 
that  offer  both  structural  load  bearing  capabilities  as  well  as  low  stiffness  configurations  during 
reconfiguration.  For  this  purpose,  SMP  filled  honeycomb  cores  have  been  proposed.  [29,  30,  31, 
25,  32] 

One  of  the  main  barriers  to  implementation  of  such  a  skin  system  is  the  method  used  to  heat  the 
SMP.  Methods  of  making  SMPs  conductive  and  thus  electrically  activated  through  joule  heating 
have  been  extensively  studied.  Systems  including  carbon  fibers,  carbon  nanotubes,  carbon  black, 
exfoliated  graphite,  metallic  nanostrands,  metallic  powders,  small  wires,  and  oxide  powders  have 
all  been  proposed.  [33,  34,  35,  36,  37]  Luo  and  Mather  were  able  to  achieve  sub-2-second 
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activation  using  continuous  carbon  nanofibers  in  an  epoxy  SMP  matrix.  The  glass  transition 
temperature,  Tg,  for  the  sample,  50  °C,  was  lower  than  many  applieations  require  however  and 
the  sample  was  very  thin,  whereas  transitioning  a  thick  sample  is  more  eomplicated.  Niekel 
coated  graphite  fibers  [33,  34]  and  stainless  steel  fibers  [34]  have  been  used  to  ereate  conduetive 
polypropylene  composites  with  some  success  however  require  sophisticated  manufacturing 
techniques.  Nickel  powders  have  been  used  similarly  to  the  presented  study  by  Leng  et  ah, 
providing  joule  heating  by  aligning  the  particles  into  ehains  using  lower  viscosity  resins  and  a 
weaker  polymer.  [35,36]  Similar  to  the  presented  work,  they  aligned  metallic  particles  using 
magnets  in  a  thin  strip  of  uncured  polymer  resin.  The  resulting  nickel  ehains  greatly  redueed  the 
electrieal  resistance  of  the  sample  and  were  shown  to  retain  their  integrity  after  several  strain 
cyeles. 

3.3.1  Particle  Chain  Spacing 

The  SMP  used  in  the  current  study  is  an  epoxy  SMP  comprising  of  EPON826  and 
JeffamineD230.  This  polymer  was  chosen  for  its  known  charaeteristics  and  tailorable  properties. 
Compositions,  eure  eycle,  and  manufaeturing  procedures  are  published  by  Xie  et  al.  [38]  The 
resistive  heating  elements  embedded  in  the  SMP  are  constructed  of  Nickel  particles,  ehosen  for 
their  resistive  heating,  magnetie,  and  eorrosion  properties.  The  particles  are  approximately  5  pm 
in  diameter  spheres  supplied  by  Sigma- Aldrieh.  Material  properties  for  both  the  SMP  and 
Nickel  particles  are  shown  in  Table  2  and  Table  3. 


Table  2:  Nickel  Particle  Properties 


Material  Property 

Value 

Density  (g/m  ) 

8.908e6 

Young’s  Modulus  (Pa) 

207e9 

Electrieal  Resistivity  (Q-em) 

6.4e-6 

Magnetic  Permeability 

1240 

Particle  Diameter  (m) 

5e-6 

Thermal  Capacity  (J/m^/C) 

0.46 

Thermal  Conduetivity  (W/m/K) 

60.7 

Table  3:  Shape  Memory  Polymer  Properties 

Material  Property 

A 

H 

OTQ 

V 

H 

CTQ 

Young’s  Modulus  (Pa) 

1.3e9  1.9e7 

Thermal  Cond.  (W/m/K) 

0.206  0.234 

Thermal  Capacity  (J/m^/K) 

1.43e6  2.45e6 

Density  (g/m^) 

1.153e6 

Tg(C) 

78 

Viscosity  25C  (Pa-s) 

8 

Viscosity  75C  (Pa-s) 

0.08 

Thermal  Expansion  Coeffieient  (m/m/K) 

55e-6 

Some  of  the  advantages  of  the  ehosen  heating  scheme  are  its  inherent  flexibility  in  material 
seleetion,  ease  of  manufaeturability,  and  tailorable  properties.  Unfortunately  sinee  the  polymer 
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has  a  relatively  low  thermal  diffusivity,  either  extremely  long  times  are  required  to  heat  the 
polymer  to  transition  or  the  material  nearest  the  heat  source  will  be  significantly  above  Tg;  both 
of  which  are  undesirable.  One  of  the  principal  advantages  of  distributed  heat  sources  such  as  the 
presented  Nickel  particle  chains  is  the  excess  thermal  energy  required  to  heat  the  polymer 
quickly  is  minimized  by  minimizing  the  distance  between  heat  sources. 

Using  a  simple  ID  thermal  model  described  by  Equations  (13)  through  (24),  the  minimum 
volume  density  of  particles  needed  to  heat  the  SMP  to  Tg  in  a  given  time  can  be  calculated. 
Particle  chains  are  modeled  as  cylinders  with  effective  radii  based  on  surface  area,  r  in  Figure  16. 
If  the  chains  are  efficiently  packed  in  the  material  and  spaced  2b  apart,  then  each  chain  is 
approximately  responsible  for  heating  a  cylindrical  area  of  radius  b  ’  from  the  center  of  the  chain. 
The  SMP  is  then  modeled  as  a  hollow  cylinder  with  inner  radius  r  and  outer  radius  b  The  outer 
radius  is  assumed  to  be  adiabatic  due  to  symmetry  while  a  uniform  heat  flux,  q^,  is  applied  to  the 
inner  radius.  The  cylinder  is  assumed  to  be  at  uniform  initial  temperature  and  thermal  properties 
of  the  polymer  are  allowed  to  vary  with  temperature;  which  has  been  shown  to  have  a  significant 
impact  on  the  calculated  thermal  energy  and  time  required  to  reach  transition.  [39]  With  an 
initial  temperature  of  25  °C  and  the  temperature  at  the  interface  between  the  chain  and  polymer 
limited  to  200  °C  to  prevent  damage  to  the  polymer,  the  time  required  for  the  temperature  b  ’ 
from  the  center  of  the  chain  to  reach  Tg  is  calculated.  Supplied  power  is  distributed  over  a  cube 
of  SMP  1  cm  by  1  cm  in  area  and  5  mm  thick. 


Figure  16:  Parameters  for  ID  Thermal  Diffusion  Between  Chains 
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Figure  17  shows  the  time  required  for  the  temperature  b  ’  from  the  particle  chain  to  reach  78  °C 
with  respect  to  the  volume  fraction  of  Nickel  particles.  An  order  of  magnitude  increase  in  the 
volume  fraction  of  Nickel  particles  results  in  only  a  12%  reduction  in  heating  time. 

Interestingly,  the  low  thermal  conductivity  of  the  polymer  and  limiting  the  temperature  at  the 
interface  to  prevent  damage  results  in  the  peak  supplied  power  having  little  effect  on  the  excess 
energy  required  to  reach  transition,  thus  only  one  curve  is  visible  in  Figure  18.  Excess  energy 
being  defined  as  the  amount  of  thermal  energy  delivered  to  the  polymer  beyond  what  would  be 
required  if  the  polymer  were  uniformly  at  Tg.  In  each  case  the  temperature  of  the  polymer  near 
the  particle  chain  reaches  200  °C  before  the  temperature  at  b  ’  reaches  Tg.  This  results  in  the 
amount  of  excess  thermal  energy  required  having  a  similar  trend  to  that  of  the  time  required  for 
transition  with  approximately  13%  less  energy  being  required  using  10  vol%  Nickel  than  1  vol% 
while  the  peak  supplied  power  has  little  to  no  effect. 
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Figure  17:  Minimum  Time  to  Reach  Tg  Given  Volume  Fraction  of  Particles  and  Heating 

Power 


Figure  18:  Excess  Energy  Required  to  Reach  Tg  Given  Volume  Fraction  of  Nickel  Particles 

A  system  with  perfectly  dispersed  chains  is  of  course  unrealistic.  The  chains  formed  in  a  real 
system  are  not  evenly  spaced,  have  non-connecting  branches  and  sections,  and  may  be  more  than 
one  particle  in  thickness.  Thus,  the  above  predictions  represent  the  minimum  power  and  particle 
fraction  required  for  a  desired  heating  time. 

3.3.2  Nickel  Particle  Chain  Formation  and  Characterization 

3.3.2. 1  Particle  Alignment  Due  to  Magnetic  Fields 

To  align  the  particles  into  chains,  randomly  dispersed  particles  in  uncured  SMP  resin  are 
subjected  to  a  magnetic  field.  Neodymium  N42SH  magnets,  serviceable  up  to  150  °C,  encased 
in  aluminum  fixtures  are  suspended  above  and  below  the  sample  as  shown  in  Figure  19  and 
Figure  24.  The  magnets  are  60  mm  apart  with  the  sample  placed  in  the  center.  Each  magnet  is 
101.6  mm  wide  and  3.175  mm  thick.  The  sample  shown  in  Figure  19  is  127  mm  in  length  and 
5mm  thick.  The  polymer  resin  sample  is  in  an  open  Teflon®  mold  which  sits  in  two  “C”  shaped 
Teflon®  guide  blocks.  The  magnetic  field  surrounding  the  sample  is  modeled  using  the  freeware 
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program  Finite  Element  Method  Magnetics  (FEMM)  version  4.2.  The  magnetic  field  within  the 
sample  has  an  average  strength  of  4.77e4  A/m,  with  a  standard  deviation  of  1.9e4  A/m  (40%). 
While  there  is  significant  variation  in  magnetic  field  strength  horizontally,  particularly  near  the 
edges,  vertically  the  field  is  nearly  uniform  with  a  standard  deviation  of  1.82e3  A/m  (3.8%). 
Placing  the  magnets  closer  together  increases  the  vertical  standard  deviation  of  field  strength 
making  the  field  less  uniform.  A  non-uniform  field  causes  a  net  non-zero  magnetic  force 
between  the  particles  in  the  sample  and  the  magnets  which  causes  the  particles  to  collect  on  the 
top  and  bottom  surfaces  rather  than  form  chains. 

Also  shown  is  the  induced  magnetic  field  between  three  vertically  aligned  Ni  particles.  The 
magnetic  polarization  created  in  the  particles  dramatically  increases  the  field  density  between 
nearby  particles.  The  three  particles  shown  in  Figure  19  are  5  pm  and  15  pm  apart.  The 
magnetic  polarization  created  in  the  particles  by  the  magnets  causes  the  particles  to  be  attracted 
to  one  another  when  aligned  vertically;  forming  chains  along  the  magnetic  field  lines  between 
the  two  magnets.  The  entire  assembly  is  placed  in  an  oven  at  100  °C  for  3  hours  to  cure, 
ensuring  good  particle  alignment  during  the  expansion  and  shrinking  that  occurs  during  curing. 
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Figure  19:  Magnetic  Field  Intensity  Surrounding  the  Sample  (Top)  and  Between  Particles 

(Bottom) 
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Figure  20:  Magnetic  Field  Intensity  and  Direction  Near  a  Nickel  Particle 


Using  FEMM,  the  attractive  or  repulsive  force  between  a  pair  of  particles  can  be  calculated. 
Particles  aligned  vertically  attract  while  particles  in  the  same  plane  aligned  horizontally  repel. 
The  magnetically  induced  attractive  force  between  two  vertically  aligned  spherical  Nickel 
particles  is  shown  in  Figure  21  as  a  function  of  the  distance  between  the  particles.  A  similar 
negative  curve  can  be  generated  for  the  repulsive  force  between  particles  in  plane. 


Figure  21:  Magnetic  Attractive  Force  Between  Vertically  Aligned  Particles 

The  net  downward  force  on  a  particle  taking  into  account  gravity  and  buoyancy  is  4.97  pN. 
Thus,  for  any  positive  vertical  movement  the  magnetic  force  between  two  vertically  aligned 
particles  must  be  higher  than  4.97  pN.  With  the  above  magnetic  setup  this  occurs  only  for 
particles  less  than  4.2  pm  apart.  At  1  vol%,  if  the  particles  are  efficiently  dispersed  there  is  an 
average  19  pm  between  each  particle.  At  19  pm,  the  magnetic  force  between  two  vertically 
aligned  particles  is  0.3 1  pN,  far  less  than  is  required.  While  the  magnetic  field  is  sufficient  to 
vertically  align  particles  in  close  proximity,  it  is  insufficient  to  properly  attract  particles  more 
than  a  diameter  apart. 

It  is  possible  to  increase  the  strength  of  the  magnetic  field  and  thus  the  induced  magnetic  field 
within  each  particle  resulting  in  stronger  magnetic  forces  between  particles.  The  distance 
between  the  magnets  pictured  in  Figure  19  can  be  decreased  or  the  strength  of  the  magnets  can 
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be  increased.  Either  of  these  methods  however  would  increase  the  magnetic  field  gradient 
experienced  by  the  sample.  A  strong  gradient  results  in  a  net  non-zero  magnetic  force  between 
the  particles  and  magnets  attracting  the  particles  to  either  the  top  or  bottom  surface  of  the 
sample.  While  not  directly  studied,  increased  particle  densities  near  the  surface  of  the  sample 
and  significant  surface  texturing  were  seen  for  magnet  spacings  less  than  40  mm  apart. 

3.3.2.2  Particle  Alignment  Due  to  Linear  Vibrations 

To  supplement  the  force  on  the  particles  generated  by  the  magnetic  field,  linear  polarized 
vibration  is  introduced  to  the  sample  to  increase  particle  mobility.  Particle  motion  in  a  vibrating 
fluid  has  been  extensively  studied  and  modelled.  [40]  Figure  22  shows  the  direction  of  forces 
generated  on  a  particle  for  any  pair  of  particles  aligned  along  one  of  the  three  cardinal  axes  when 
subjected  to  a  polarized  vibrating  fluid. 


Figure  22:  Gravity,  Buoyancy,  Vibration,  and  Magnetic  Forces  on  a  Suspended  Particle 

In  Figure  22  Fb  is  the  buoyancy  force,  Fg  is  the  force  due  to  gravity,  and  V  and  M  are  forces 
induced  by  vibration  and  the  magnetic  field  respectively  with  respect  to  the  nearby  particle  noted 
in  the  subscript.  Tangential  vibratory  forces  such  as  that  created  between  particles  Pq  and 
Picause  neighboring  particles  to  repel  while  normal  vibrations  such  as  that  between  Pq  and  P2 
cause  particles  to  attract.  If  the  central  particle  Pq  were  surrounded  by  particles  on  the  negative 
axes  as  well  as  the  particles  shown  in  Figure  22  (six  particles  total)  and  each  were  the  same  size, 
shape,  and  distance  from  particle  Pg,  the  resulting  forces  due  to  the  vibrating  fluid  and  magnetic 
field  would  cancel  and  the  net  force  on  the  particle  would  simply  be  the  difference  between  the 
force  due  to  gravity  and  the  buoyancy  force.  A  perfectly  uniform  distribution  of  particles  would 
then  result  in  the  particles  simply  settling.  Non-uniform  particle  shape  and  size  and  varying 
distances  between  particles  is  believed  to  result  in  a  net  non-zero  force  on  each  particle  causing 
them  to  preferentially  align  along  magnetic  field  lines.  Since  all  forces  in-plane  cause  particles 
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to  repel,  the  vibration  frequency  and  amplitude  and  strength  of  the  magnetic  field  can  be 
optimized  such  that  the  distribution  of  chains  is  relatively  uniform. 

Lyubimov  et  al.  derived  the  below  equations  for  the  force  experienced  by  a  pair  of  particles  in  a 
vibrating  fluid  dependent  upon  their  orientation  with  respect  to  the  vibration. 
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Where  a  and  cu  are  the  vibration  amplitude  and  frequency,  r  is  the  particle  radius,  pp  is  the 
particle  density,  p/is  the  fluid  density,  Vp  is  the  particle  volume,  8  is  the  distance  between  the  two 

particles,  ^  is  a  unit  normal  vector  connecting  the  two  particles,  ^  is  the  gradient  of  the  mean 

squared  pulsational  velocity  averaged  over  the  vibration  period,  2  *  ^  -  cos  6  ^  p  •  ^ 
force  experienced  by  each  particle.  Particles  normal  to  the  flow  (0  =  7r/2)  experience  an 
attractive  force  while  particles  aligned  parallel  to  the  flow  (0  =  0)  repel. 

As  predicted  by  Equations  (25)  and  (26),  particles  normal  to  the  flow  experience  a  positive  force 
causing  them  to  attract  while  the  magnetic  forces  and  forces  due  to  parallel  flow  cause  particles 
to  repel.  This  interplay  between  attractive  and  repulsive  forces  in  plane  causes  banding  to  occur 
in  the  sample  where  parallel  rows  of  particles  accumulate  perpendicular  to  the  direction  of 
vibration.  This  alignment  of  particle  chains  has  been  seen  experimentally  by  Lyubimov  for  2- 
degree-of- freedom  systems  in  which  a  single  layer  of  particles  are  subjected  to  vibrations  in  a 
fluid.  The  effect  is  also  observable  in  samples  from  the  presented  work  many  particles  thick  as 
shown  in  Figure  23. 


Figure  23:  Backlit  10-mm  Thick  Sample 


27 

Approved  for  public  release;  distribution  unlimited. 


In  Figure  23,  particles  are  dark  with  a  light  background.  Distinct  banding  spaced  approximately 
1  mm  apart  is  highlighted  with  red  arrows.  The  sample  was  subjected  to  both  the  magnetic  field 
pictured  in  Figure  19  and  50  pm  linear  vibrations  at  300  Hz.  Slight  particle  alignment 
perpendicular  to  the  direction  of  vibration  in  a  1 -mm- thick  sample  subjected  only  to  vibration 
(no  magnetic  field  was  applied)  was  also  observed.  While  this  and  Figure  23  indicate  that 
Equations  (25)  and  (26)  are  applicable  to  the  presented  case,  they  are  as  of  yet  unverified  for 
vibrational  amplitudes  larger  than  the  diameter  of  the  particles  being  modelled  and  thus  are 
provided  here  as  a  qualitative  reference.  It  should  also  be  noted  that  the  50  pm  300  Hz  vibration 
is  applied  to  the  mold.  No  attempt  to  measure  the  motion  of  the  fluid  was  performed.  It  is  also 
believed  that  fluid  motion  near  the  bottom,  sides,  and  open  top  surface  varies  greatly.  Imperfect 
particle  shape  and  size  differences  are  also  factors  in  the  resulting  particle  dynamics.  The 
complexity  of  factors  influencing  particle  motion  is  such  that  modeling  the  system  dynamics 
directly  is  beyond  the  scope  of  the  presented  work.  Thus,  the  presented  theories  are 
predominantly  qualitative. 

3.3.3  Sample  Fabrication 

The  experimental  setup  used  to  achieve  particle  alignment  is  shown  in  Figure  24.  The  300  Hz 
shaker  motor  signal  is  generated  by  a  function  generator.  The  fulcrum  connecting  the  shaker  and 
the  mold  reduces  the  stroke  of  the  shaker  10:1  resulting  in  the  sample  experiencing  an  amplitude 
of  approximately  50  pm.  After  curing,  samples  are  lightly  sanded  and  copper  mesh  electrodes 
are  applied  to  the  top  and  bottom  surfaces  using  the  same  epoxy  SMP  with  10  vol%  nickel 
particles.  The  electrodes  are  clamped  between  two  thin  sheets  of  Teflon  and  aluminum  plates 
during  curing  ensuring  good  contact. 


Figure  24:  Experimental  Magnetic  and  Vibration  Setup 
3.3.4  Experimental  Results 
3.3.4.1  Particle  Alignment 

Using  the  above  method,  nickel  particle  chains  can  be  readily  formed  in  epoxy  resin  during 
curing.  The  magnetic  and  vibratory  forces  are  sufficient  to  generate  particle  motion  such  that  no 
delay  is  required  between  mixing  and  curing  the  sample.  The  particles  align  in  the  oven  during 
the  first  few  minutes  of  the  cure  cycle  before  any  significant  polymerization.  Figure  25  shows 
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the  alignment  achieved  in  a  typical  sample  with  1  vol%  Nickel  particles.  While  there  are  several 
chains  that  terminate  in  the  middle  of  the  sample  and  the  chains  vary  in  thickness,  there  is 
generally  good  alignment  and  distribution. 


Figure  25:  lOX  (Left)  and  SOX  (Right)  View  of  1  Vol%  Nickel  Particle  Chains 
3.3.4.2  Electrical  Properties 

As  a  measure  of  the  quality  of  particle  alignment  and  the  heating  capability  of  the  nickel  chains, 
the  resistance  of  samples  at  various  temperatures  was  measured  and  the  resistivity  of  the  material 
calculated.  Figure  26  shows  the  results  for  a  1  vol%  Nickel  particle/epoxy  SMP  composite.  At 
ambient  temperature  the  average  resistivity  is  57  Q-mm  with  very  little  variation  between 
samples.  As  the  temperature  increases  the  resistivity  of  the  samples  increases  until  unreadably 
high  resistance  values  are  attained  above  160  °C. 


Figure  26:  Average  Resistivity  and  Standard  Deviation  of  12  Samples  Allowed  to 

Thermally  Expand 
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The  more  than  seven  order  of  magnitude  inerease  in  resistivity  ean  possibly  be  attributed  to 
several  faetors,  the  dominant  being  thermal  expansion.  Between  ambient  and  160  °C  the 
samples  expand  by  approximately  0.8%.  It  is  postulated  that  this  is  sufficient  for  the  majority  of 
particles  to  be  pulled  away  from  one  another,  ceasing  to  physically  touch,  resulting  in  an  open 
circuit.  Indeed,  if  a  sample  is  placed  in  a  vise  in  an  over  at  120  °C  and  gradually  subjected  to 
compressive  strains  the  resistivity  of  the  samples  reduce  to  ambient  temperature  levels  as  shown 
in  Figure  27. 
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Figure  27:  Average  Resistivity  of  Samples  at  120  °C  Subjected  to  Compressive  Strain 

Measured  resistivity  of  samples  of  varying  nickel  contents  are  shown  in  Figure  28.  Samples  with 
0.75  vol%  nickel  and  below  have  significantly  higher  resistivity  at  room  temperature  than  those 
with  1  vol%  or  higher.  While  nickel  content  is  believed  to  have  a  significant  effect  on  sample 
resistivity,  the  cause  of  the  difference  in  initial  resistivity  above  and  below  1  vol%  is  unclear. 

All  volume  fractions  however  increase  in  resistivity  with  increasing  temperature,  previously 
shown  to  be  the  result  of  thermal  expansion. 
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Figure  28:  Resistivity  Dependence  on  Temperature  and  Vol%  Nickel  Particles 
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3.3.4.3  Thermal  Distribution 

Regardless  of  the  eause,  the  inerease  in  resistivity  results  in  two  advantages.  Due  to  the  variation 
in  the  quality  and  distribution  of  the  Nickel  particle  chains  the  samples  are  not  perfectly  uniform 
in  their  electrical  properties.  There  exist  pathways  that  are  more  conductive  than  others  through 
the  sample,  causing  current  to  preferentially  flow  through  certain  chains  and  not  others.  If 
unaltered,  this  would  cause  the  samples  to  heat  up  non-uniformly  drastically  increasing  the  time 
and  energy  required  to  reach  transition.  Since  the  resistance  of  the  sample  increases  with 
temperature  (due  to  local  thermal  expansion),  areas  of  high  temperature  have  higher  resistance 
than  those  of  lower  temperature.  As  current  flows  through  a  certain  chain  or  set  of  chains,  the 
temperature  increases.  As  the  resistance  increases  with  increasing  temperature;  the  path  of  least 
resistance  shifts  to  another  area.  As  current  begins  to  flow  through  this  new  area;  the 
temperature  increases  and  the  cycle  continues.  This  results  in  the  sample  heating  more  evenly 
and  quickly  and  can  clearly  be  seen  in  the  series  of  infra-red  images  in  Figure  29.  The  top  image 
shows  the  sample  just  after  current  has  been  applied  where  the  left  most  particle  chain  is  slightly 
warmer  than  the  surrounding  sample.  After  0.55s,  the  temperature  of  the  chain  is  high  enough 
that  the  path  of  least  resistance  has  moved  to  the  right,  Figure  29  middle.  After  0.935s  the  SMP 
near  the  particle  chains  on  the  left  side  of  the  sample  has  reached  Tg  and  the  sample  in  this  area 
has  become  non-conductive.  Figure  30  shows  a  magnified  section  of  the  sample  shown  in  Figure 
29  at  0.935s  overlaid  with  an  optical  image  of  the  sample.  Areas  of  increased  temperature 
clearly  align  with  areas  with  high  densities  of  particle  chains. 
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Figure  29:  Infrared  Images  of  30V,  O.OIA  Applied  to  8mm  Tall  by  2mm  Thick  Sample. 
0.055s  (Top),  0.55s  (Middle),  and  0.935s  (Bottom)  (Temperature  scale  in  °C) 


70 

60 

50 


Figure  30:  Infrared  (Top),  Optical  (Bottom),  and  Overlay  (Middle)  Images  of  Nickel  Chain 

Heating  (Temperature  scale  in  °C) 

The  second  advantage  of  the  composite  increasing  in  resistivity  with  increased  temperature  is  the 
ability  to  use  the  resistance  of  the  sample  as  a  temperature  sensor.  Rather  than  requiring  a 
separate  thermocouple  to  determine  when  the  desired  temperature  has  been  reached,  the 
resistance  of  the  sample  can  be  used  to  estimate  temperature. 
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3.3.5  Nickel  Chain  Heating  Elements  Conclusions 

A  method  for  aligning  magnetic  particles  in  fluids  with  high  viscosities  has  been  demonstrated 
leading  to  an  effective  technique  for  quickly  and  uniformly  heating  shape  memory  polymer.  It 
has  been  shown  that  sub-two  second  activation  times  are  achievable  with  very  low  volume 
fractions  (<1.0%)  of  conductive  particles  minimizing  the  added  weight  associated  with  many 
heating  methods.  The  low  volume  fractions  are  made  possible  by  aligning  nickel  particles  into 
chains  under  a  magnetic  field  and  polarized  vibration.  The  requirement  that  the  particles  be 
magnetic  however  limits  the  selection  of  particle  materials  and  precludes  many  traditional 
dopants  such  as  carbon  black.  Such  dopants  could  however  be  used  to  increase  the  diffusivity  of 
the  composite,  further  reducing  the  transition  time  by  decreasing  temperature  gradients  in  the 
polymer.  The  composite’s  increase  in  resistivity  with  increasing  temperature  was  characterized 
and  utilized  experimentally  to  aid  in  uniformly  heating  the  sample  and  as  a  thermoelectric  switch 
signaling  when  the  sample  had  reached  the  transition  temperature.  Significant  differences  in  the 
time  and  power  required  to  reach  transition  were  measured  and  are  believed  to  be  due  to 
variations  in  the  quality  of  particle  alignment  within  the  samples.  With  optimization  of 
manufacturing  techniques  the  excess  energy  required  over  the  thermal  capacity  of  the  polymer 
should  be  able  to  be  reduced,  making  the  presented  heating  method  even  more  attractive  for 
applications.  The  presented  SMP  heating  method  adds  minimal  additional  weight,  is  simple  to 
manufacture,  results  in  uncomplicated  designs,  is  low  cost,  and  results  in  a  relatively  uniform 
temperature  making  it  ideal  for  many  shape  memory  polymer  applications. 

3.4  Material  Characterization 

Characterization  of  the  material  used  in  the  presented  study  is  detailed  in  the  published  interim 
report:  FA8650-09-D-3944  Interim  Report.  Below  is  a  brief  synopsis  of  the  conclusions. 

3.4.1  SMP  Characterization 

The  SMP  infill  used  to  experimentally  validate  the  model  is  a  two  part  epoxy  consisting  of  7.28  g 
EPON  826  to  2.3  g  Jeffamine  D-230  published  by  Xie  and  Rousseau  [41].  The  hand-mixed  resin 
is  placed  in  a  medium  vacuum  to  remove  air  bubbles  and  cured  in  an  oven  at  100  °C  for  3  hours 
in  an  open  Teflon®  mold.  Dogbone  samples  were  cut  from  sheets  of  neat  SMP  according  to 
ASTM  standard  D638  and  tested  in  a  screw-driven  MTS  ATest/lL  tensile  machine  with  an  MTS 
50  Ibf  load  cell.  An  MTS  634-1  lE-24  laser  extensometer  was  used  to  measure  strain  with  two 
thermocouples  measuring  air  and  sample  temperature.  Heat  was  supplied  to  a  closed  insulated 
chamber  via  a  variable  voltage  source,  electric  resistance  in-line  heater,  and  compressed  house 
air.  The  sample  was  allowed  to  reach  equilibrium  at  each  temperature  for  approximately  five 
minutes  before  testing.  Samples  6mm  wide  by  2.5-mm  thick  with  50-mm  gage  lengths  were 
tested  at  a  strain  rate  of  approximately  1.5  "''Vmm-  The  modulus  of  the  polymer  was  then 
calculated  at  various  temperatures  and  fit  with  a  modified  Heaviside  function  to  determine  the 
cold  and  hot  state  moduli  and  glass  transition  temperature,  Tg.  Results  are  listed  in  Table  4. 
Experimental  error  values  where  appropriate  to  report  represent  95%  confidence  intervals.  Shear 
tests  on  neat  SMP  were  also  conducted  using  the  fixture  and  method  described  below. 
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Table  4:  Epoxy  SMP  Experimental  Results 


Material  Property 

Young’s  Modulus  <  Tg 
Transition  Temperature 
Young’s  Modulus  >  Tg 
Shear  Modulus  <  Tg 
Shear  Modulus  >  Tg 


Value 


1309±120  MPa 
78±4  °C 
19.04±3.1  MPa 
1270  MPa 
1.06  MPa 


3.4.2  Honeycomb  Characterization 


The  honeycomb  utilized  in  the  study  is  aluminum  alloy  3003  from  McMaster-Carr  with  the 
dimensions  listed  in  Table  5. 

Table  5:  Experimental  Honeycomb  Geometry 

Dimension  Value 

a  5.25  mm 

1  9.33  mm 

d  0.044  mm 

c  6.39  mm 

0  0.799  rad 

Eh  73.1  GPa 


Samples  were  cut  using  scissors  into  sheets  127  mm  wide  by  254  mm  long  at  orientations  of  a  = 
0,  Vi2,  %,  %,  and  "/2  and  tested  using  the  same  MTS  load  frame  and  load  cell  as  the  SMP 
samples  at  a  strain  rate  of  1.8  A  laser  extensometer  was  used  to  measure  strain  with  a  gage 
length  of  approximately  50  mm.  Resulting  values  of  Young’s  modulus  at  various  angles  are 
listed  in  Table  6  as  well  as  their  associated  95%  confidence  intervals. 

Table  6:  Experimental  Young’s  Moduli  for  Empty  Honeycomb 


Angle  (a) 

Young’s  Modulus  (Pa) 

0° 

5.55E4±  1.17E4 

15° 

4.72E5  ±  0.58E5 

30° 

7.83E5±  1.30E5 

45° 

6.37E5  ±  0.72E5 

SO 

o 

o 

2.01E4±  1.83E4 

Several  models  available  in  the  literature  capable  of  predicting  the  in-plane  moduli  of 
honeycomb  are  listed  in  Table  10.  Each  model  is  used  to  predict  the  tensile  modulus  of  an 
Aluminum  Alloy  3003  (Young’s  modulus  73.1  GPa,  shear  modulus  25  GPa)  honeycomb  with 
the  dimensions  listed  in  Table  5.  Since  the  chosen  shape  is  not  a  standard  honeycomb  where  a  is 
equal  to  /  and  0  equal  to  "/s;  the  moduli  in  the  a=0  and  %  directions  are  not  equal,  as  shown  in 
Table  7. 
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Table  7:  Comparison  Between  Model  Predictions  of  Young’s  Modulus  for  Empty 

Honeycomb 

Young’s  Modulus  (kPa) 

_ Eo _ E,i/2 _ 


Beblo 

26.24 

8.945 

Masters  and  Evans 

26.24 

8.945 

Abd  El-Sayed,  et  al 

26.24 

8.945 

Gibson  and  Ashby  (simple) 

26.24 

8.946 

Gibson  and  Ashby  (full) 

209.8 

71.54 

Experimental  Values 

55.5 

20.1 

With  the  exception  of  the  unabridged  model  by  Gibson  and  Ashby,  all  of  the  models  give  nearly 
identical  results  underestimating  the  experimentally  determined  values. 

3.4.3  Composite  Characterization 

3.4.3.1  Experimental  Composite  Characterization 

Honeycomb  samples  were  cut  and  etched  in  a  23.5%  nitric  acid  bath  for  10  minutes  to  clean  the 
honeycomb  and  promote  adhesion  between  the  aluminum  and  the  polymer.  Samples  were  then 
rinsed  and  dried  in  an  oven.  Epoxy  SMP  was  prepared  as  above  and  a  thin  coating  applied  to  the 
bottom  of  a  Teflon  mold.  The  aluminum  honeycomb  was  then  placed  in  the  mold  and  the 
assembly  cured  in  an  oven  at  100°C  for  an  hour.  Every  other  row  of  cells  was  then  filled  with 
SMP  and  the  assembly  cured  for  another  hour.  Finally,  the  remaining  cells  were  fdled  with  SMP 
and  the  completed  mold  cured  at  100°C  for  an  additional  three  hours.  This  process  was  found  to 
reduce  damage  due  to  shrinkage  during  curing. 

Young’s  modulus  was  measured  in  the  a=0  and  %  directions  using  127  mm  by  254  mm 
specimens,  an  MTS  2.2  kip  screw  driven  load  frame  with  a  2,250  Ibf  load  cell,  and  a  laser 
extensometer  at  1 .27  ”™/min  strain  rate.  Above  Tg  tests  were  conducted  using  a  custom  insulated 
thermal  camber,  in-line  electrical  resistance  heater,  house  compressed  air,  and  two 
thermocouples  measuring  air  and  sample  temperature. 

Table  8  shows  experimental  Young’s  moduli  predictions  for  filled  honeycomb  above  and  below 
Tg.  The  range  of  values  of  stiffness  in  various  directions  is  significantly  reduced  compared  to 
empty  honeycomb  as  the  addition  of  SMP  has  an  averaging  effect  on  the  stiffness  of  the 
composite.  The  difference  between  the  analytically  predicted  Young’s  modulus  in  the  a=0  and 
"/i  directions  and  the  experimentally  determined  values  are  less  than  6.5%  and  4.3% 
respectively. 
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Table  8:  Experimental  Composite  Young’s  Modulus 


Direction  (a) 


(Mp‘a) 


(mA) 


0 

0 

"/2 

0+"/2 


33.9±4.1  2190±110 

NA  NA 

11.8±4.2  204O±140 

NA  NA 


Shear  specimens  were  100  mm  square  with  holes  drilled  %  inch  apart  in  a  pattern  along  all  four 
sides  such  that  the  aluminum  honeycomb  was  not  damaged.  Specimens  were  loosely  bolted  to  a 
picture  frame  style  fixture  while  shear  deformation  was  measured  using  digital  image 
correlation.  Tests  were  conducted  on  an  MTS  2.2  kip  screw  driven  load  frame  with  a  2,250  Ibf 
load  cell.  Above  Tg  tests  were  conducted  with  a  custom  insulated  chamber,  in-line  electric 
heater,  house  compressed  air,  and  two  thermocouples  measuring  air  and  sample  temperature.  A 
typical  picture  frame  fixture  mounted  as  a  diamond  in  the  load  frame  requires  that  the  shear 
modulus  of  the  material  in  the  XY  and  YX  directions  be  equal.  Since  honeycomb  is  anisotropic 
in  shear,  tests  were  conducted  in  a  custom  shear  test  fixture  simulating  pure  shear.  The  left  side 
of  the  fixture  at  both  the  top  and  bottom  were  pinned  to  the  load  frame  eliminating  translation  in 
the  vertical  and  horizontal  directions  while  allowing  rotation.  The  bottom  right  pin  of  the  fixture 
was  free  while  the  top  right  pin  was  attached  to  the  crosshead  of  the  load  frame.  This  setup 
allows  characterization  of  both  the  XY  and  YX  shear  properties  of  the  material. 

Table  9:  Experimental  Composite  Shear  Modulus 


>Tg 

<Tg 

Direction  (a) 

(MPa) 

(MPa) 

o,% 

13.9±3.5 

1190±130 

%,o 

13.0±2.7 

1130±13O 

3.4.3.2  FEA  Composite  Characterization 

Extensive  FEA  modeling  of  a  unit  cell  of  filled  honeycomb  was  performed,  the  details  of  which 
are  presented  in  the  interim  report  with  a  brief  synopsis  below.  Periodic  boundary  conditions 
were  applied  to  a  unit  cell  of  the  composite  in  Abaqus  using  the  materials  and  material  properties 
listed  previously,  the  results  of  which  are  used  to  investigate  the  mechanics  of  the  honeycomb 
cell  and  validate  the  experimental  and  analytical  results  and  assumptions.  Test  cases  modeling 
the  SMP  infill  in  both  the  hard  and  soft  state  were  conducted  in  the  X,  Y,  XY,  and  YX  directions 
with  periodic  boundary  conditions.  Figure  3 1  shows  stress  in  the  0°  direction  when  the  sample  is 
subjected  to  strain  in  the  0°  direction.  The  stress  in  the  infill  along  the  slanted  members  in 
Figure  3 1  ranges  from  -0. 12  to  0.2  Pa  nearly  linearly  with  the  stress  near  the  center  of  the  beam 
approaching  zero.  Such  distribution  validates  the  same  assumption  employed  in  the 
development  of  the  analytical  model  presented  subsequently. 
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Figure  31:  Stress  in  the  0°  Direction  With  Soft  Infill  (Pa) 


Figure  32  shows  the  stress  in  the  infill  in  the  nil  direction  when  subjected  to  strains  in  the  0° 
direction.  The  stress  along  the  slanted  members  is  entirely  negative  and  nearly  uniform, 
correlating  well  with  the  model  assumptions  shown  in  Figure  52. 
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Figure  32:  Stress  in  the  nil  Direction  With  Soft  Infill  (Pa) 

The  stress  field  as  predicted  by  FEA  can  also  be  compared  to  experimental  digital  image 
correlation  (DIG)  results  as  shown  in  Figure  33.  For  clarity  an  orange  hexagon  is  overlaid  on  the 
DIG  results  showing  the  approximate  location  of  the  honeycomb.  Areas  of  high  and  low  stress 
between  the  two  results  correlate  well,  validating  the  homogenization  scheme  and  periodic 
boundary  conditions  used  in  the  FEA. 
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Figure  33:  Top:  FEA  Von-Mises  Stress,  Bottom:  DIC  Experimental  Von-Mises  Stress 
3.5  System  Characterization 
3.5.1  Fabrication  Procedure 

Fabrication  of  the  composite  panels  is  a  multistep  process  ensuring  electrical  connectivity  to 
each  cell  while  isolating  individual  cells  from  neighboring  cells  both  electrically  and  thermally. 
Each  step  in  the  process  is  detailed  below. 

Step  1 

Aluminum  3003  honeycomb  purchased  from  McMaster-Carr  with  dimensions  listed  in  Table  5  is 
cut  to  the  desired  size  using  scissors.  The  samples  are  then  etched  for  10  minutes  in  a  Nitric 
Acid  solution:  1  part  Fisher  Scientific  Nitric  Acid  (Certified  ACS  Plus  70%  HN03)  to  3  parts  tap 
water  (pH  8.0).  Samples  are  then  thoroughly  rinsed  and  allowed  to  air  dry.  This  process  helps  to 
clean  the  honeycomb  of  residue  and  oxidation. 
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Figure  34:  Nitric  Acid  Etching  of  Aluminum  Honeycomb 


Step  2 

To  electrically  isolate  the  honeycomb  two  coats  of  neat  Epoxy  SMP  are  deposited  on  all  surfaces 
of  the  honeycomb  using  a  cotton  tipped  applicator.  The  sample  is  placed  in  a  chemical  hood 
overnight  at  room  temperature  and  cured  in  a  100°C  oven  for  three  hours  the  following  morning. 
This  allows  the  SMP  to  partially  set  preventing  excessive  loss  during  the  initial  stages  of  curing. 
Alternatively,  the  SMP  could  be  applied  using  an  air  bmsh  or  the  sample  dipped  in  a  large 
container. 


Figure  35:  Aluminum  Honeycomb  Coated  with  SMP  (Left)  and  the  Resulting  Cellular 
Separation  (Right)  Visible  via  a  Backlit  Sample 


Step  3 

Shrinkage  of  the  SMP  during  curing  can  cause  delamination  between  the  SMP  and  the 
honeycomb  with  large  samples.  To  prevent  this,  the  honeycomb  is  filled  in  stages.  Initially  a 
thin  layer  of  the  desired  SMP/Nickel  composite  is  spread  on  the  bottom  of  the  mold  and  cured  in 
a  100°C  oven  for  1.5  hours  gluing  the  honeycomb  to  the  mold.  This  allows  every  other  row  of 
the  honeycomb  to  be  filled  and  cured  (with  or  without  a  magnetic  field  and/or  vibration) 
allowing  for  shrinkage  and  preventing  delamination. 
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Figure  36:  Honeycomb  Glued  to  Mold  Sealing  the  Bottom  of  Each  Cell  (Left)  and  a  Filled 

Sample  (Right) 


Step  4 

To  electrically  isolate  each  cell,  channels  6.35  mm  wide  and  1.5  mm  deep  are  milled  along  the 
top  and  bottom  surfaces  of  the  sample  following  the  honeycomb  providing  electrode  attachment 
points.  These  channels  are  then  fdled  with  neat  SMP. 


Figure  37:  Electrode  Attachment  Points  Before  (Left)  and  After  (Right)  Being  Backfilled 

with  SMP 


Step  5 

Sheets  of  100  by  100  copper  mesh  from  McMaster-Carr  are  applied  to  the  top  and  bottom 
surfaces  of  the  sample  using  10  vol%  Nickel/SMP.  The  sample  is  clamped  between  two  metal 
plates  of  similar  size  to  the  sample  protected  by  thin  Teflon  sheets  to  prevent  the  sample  from 
bonding  to  the  plates. 
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Figure  38:  Sample  with  Copper  Mesh  Adhered  to  Both  Surfaces 


Step  6 

Individual  electrodes  are  then  fabricated  by  milling  channels  XX  mm  wide  by  YY  mm  deep  on 
the  top  and  bottom  surfaces  of  the  sample  through  the  copper  mesh  following  the  honeycomb.  A 
vertical  column  pattern  is  traced  on  one  side  and  a  horizontal  row  pattern  is  traced  on  the 
opposite  side.  This  allows  any  single  cell  to  be  addressed  using  the  standard  row/column 
method. 


Figure  39:  Milling  of  Row  (Left)  and  Column  (Right)  Electrodes 


Step  7 

Finally,  the  traces  delineating  the  electrodes  are  fdled  with  neat  SMP  any  required  holes  or  other 
modifications  to  the  sample  are  made  for  assembly  with  the  test  fixture.  Electrical  leads  are 
soldered  to  each  row  and  column  of  copper  mesh  connecting  the  sample  to  a  voltage  source  and 
the  sample  is  painted  with  a  speckle  pattern  for  DIC.  The  final  sample  has  48  columns  and  26 
rows. 
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Figure  40:  Final  Sample  with  Electrode  Leads  Before  (Left)  and  After  (Right)  a  DIC 

Speckle  Pattern  is  Applied 


3.5.2  System  Thermal  Characterization 


To  demonstrate  the  functionality  of  the  system,  samples  similar  to  those  pictured  in  Figure  71  are 
connected  to  several  National  Instruments  NI-948 1  switch  modules  controlled  by  LabView®  and 
a  voltage  source,  see  Figure  41.  It  is  desired  to  supply  samples  with  as  much  power  as  possible 
resulting  in  the  quickest  possible  heating  times.  To  prevent  damage  to  the  modules,  voltage  and 
current  are  limited  to  41.8  V  and  3.1  A.  At  low  resistance  values,  current  limits  the  amount  of 
power  that  can  be  supplied  to  the  sample.  At  high  resistance  values,  voltage  limits  the  power 
that  can  be  supplied  to  the  sample.  Peak  power  is  attained  when  both  voltage  and  current  are  at 
their  maximum  allowable  values  when  the  sample  has  a  resistance  of  13.5  Q,  Figure  42.  Iso¬ 
voltage  and  current  lines  are  plotted  against  sample  resistance  with  40  V  and  3  A  highlighted. 
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Figure  41:  System  Thermal  Characterization  Test  Procedure 
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Figure  42.  Power  Supplied  to  the  Sample  Given  Current  and  Voltage  Limitations 

Experimentally,  the  voltage  and  current  can  be  continuously  varied  to  deliver  maximum  power  to 
the  sample.  When  the  voltage  and  current  reach  levels  predetermined  to  correspond  to  a 
sufficiently  high  resistance  indicating  the  sample  has  reached  Tg,  the  current  is  turned  off.  For 
the  example  case  shown  in  Figure  43,  current  is  applied  to  sequential  rows  of  cells  of  a  10-  by 
10-cell  panel  similar  to  pattern  four  in  Figure  71.  Initially,  maximum  power  is  achieved  by 
maximizing  the  current  with  a  low  voltage.  As  the  temperature  and  resistance  of  the  sample 
increase,  the  voltage  increases  until  the  limit  of  the  voltage  source  is  reached  at  37  seconds  at 
which  point  the  current  begins  to  decrease.  When  the  current  reaches  0.4  A  at  84  seconds,  power 
is  transferred  to  the  row  two  rows  above  the  first. 
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Figure  43.  Demonstration  of  Using  Increase  in  Resistance  to  Estimate  Temperature 

The  total  electrical  power  delivered  to  first  row  is  2125  J  while  1523  J  is  delivered  to  the  second 
row.  The  difference  is  likely  attributable  to  differences  in  the  quality  of  particle  alignment  and 
quality  of  the  electrodes.  Theoretically,  increasing  the  temperature  from  ambient  to  80  °C  of 
samples  16.4  cm^  and  6.6  mm  thick  should  require  only  907  J.  From  previous  calculations,  if 
20%  (181  J)  of  the  excess  energy  can  be  attributed  to  the  volume  fraction  of  particles.  Figure  18, 
the  remaining  1037  J  (49%)  for  the  first  row  and  435  J  (29%)  for  the  second  row  is  due  to  less 
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than  ideal  particle  chains.  It  should  be  noted  that  the  10.8  cm  rows  used  here  are  significantly 
larger  than  the  0.5  cm  samples  used  in  the  previous  analysis. 

3.5.3  System  Pressure  Characterization 

To  verify  the  heating  pattern  models  developed  in  Abaqus  and  the  potential  utility  of  the  system, 
a  full  scale  panel  was  constructed  meeting  MAS  program  dimensions.  The  panel  was 
constructed  using  the  procedure  outlined  previously  where  each  cell  is  individually  addressable. 
Various  heating  patterns  were  then  subjected  to  various  pressures  using  the  custom  manufactured 
pressure  chamber  in  Figure  44. 


A  Pressure  Chamber 


SCAlEia  WBGHT:  SHEET!  Of  1 


Figure  44:  15”  by  20”  Mold  (Left)  and  Pressure  Chamber  (Right)  Schematics 

A  full  size,  fully  instrumented  panel  was  subjected  to  an  increasing  uniform  pressure  load  with 
various  heating  patterns.  Full  field  temperature  and  deflection  were  measured  and  with  an  IR 
camera  and  DIC  respectively.  It  should  be  noted  that  while  the  below  results  are  presented  as 
quasi-static,  the  temperature  distribution  of  the  sample  is  constantly  evolving  due  to  thermal 
conductance,  as  shown  in  Figure  45.  Areas  of  cells  at  significantly  lower  temperature  are  also 
present.  While  each  cell  is  individually  conductive,  due  to  the  row/column  method  in  fabricating 
the  electrodes  to  reduce  complexity  (74  leads  vs.  1249  leads)  certain  cells  with  higher  resistances 
receive  little  current  resulting  in  little  temperature  change. 
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Figure  45:  Infrared  Images  of  an  8  Column  Heating  Pattern  at  the  Beginning  (Left)  and 

End  (Right)  of  a  Pressure  Test 

For  each  heating  pattern,  current  was  applied  to  the  desired  cells  raising  their  temperature  before 
pressure  was  applied.  The  resulting  thermal  strains  with  no  pressure  behind  the  sample  are 
shown  in  Figure  46  for  an  8  column  heating  pattern.  Areas  of  high  strain  clearly  correlate  well 
with  areas  of  high  temperature  shown  in  Figure  45. 
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Figure  46:  Thermal  Strain  on  the  Sample  Before  Pressure  is  Applied 

Each  heating  pattern  was  subjected  to  pressures  from  0  to  3.0  psi  in  0.5  psi  increments, 
exceeding  the  2.78  psi  threshold  of  the  MAS  program.  As  expected,  areas  of  higher  temperature, 
thus  lower  stiffness,  bear  higher  strain  under  pressure  as  shown  for  the  8  column  heating  pattern 
in  Figure  47. 
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Figure  47:  Major  Strain  on  the  Sample  at  3.0  psi 


Finally,  the  out  of  plane  displacement  of  the  sample  was  measured,  Figure  48,  with  the  center 
point  deflection  of  the  0,  4,  and  8  column  heating  patterns  plotted  with  respect  to  pressure  in 
Figure  49. 


Figure  48:  Z  Displacement  of  the  Sample  at  3.0  psi 
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Figure  49:  Center  Point  Deflection  of  a  Full  Size  Panel  Under  Uniform  Pressure  with 

Various  Columnar  Heating  Patterns 

As  detailed  in  Section  3.7.2,  FEA  analysis  for  various  heating  patterns  can  be  performed  yielding 
effective  in  plane  stiffnesses  that  can  be  used  to  model  the  behavior  of  the  skin  in  various 
environments.  FEA  of  an  8  column  heating  pattern  shown  in  Figure  45  results  in  Young’s 
moduli  of  2.63E8  and  1.02E9  Pa  in  the  X  and  Y  directions  respectively.  Flat  plate  theory  then 
suggests  the  center  point  deflection  of  the  sample  under  uniform  pressure  should  be 
approximately  double  the  observed  experimental  deflections  in  Figure  49.  While  there  are  many 
contributing  factors  such  as  sample  fabrication,  it  is  likely  the  areas  of  lower  than  desired 
temperatures,  and  thus  higher  modulus,  result  in  a  stiffer  than  desired  sample  and  lower 
deflections.  Although  only  columnar  patterns  are  presented  for  brevity,  virtually  unlimited 
thermal  patterns  are  attainable  with  the  described  system,  such  as  the  rows  pictured  in  Figure  50 
or  the  patterns  analyzed  performing  system  optimization  shown  in  Figure  71. 
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Figure  50:  Infrared  Images  of  a  4  Row  Heating  Pattern  Temperature  Profile 

Finally,  it  should  also  be  noted  that  the  presented  experimental  investigation  does  not  and  was 
not  intended  to  meet  MAS  requirements.  The  experimental  panel  was  fabricated  with  available 
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techniques  and  honeycomb  geometry.  The  test  successfully  demonstrates  the  concept  of  cellular 
heating  and  directionally  controllable  stiffness.  A  geometrically  optimized  solution  meeting  all 
MAS  requirements  is  presented  subsequently. 

3.6  Honeycomb  Modeling  and  Optimization 

3.6.1  Modeling  Introduction 

Since  analytical  models  describing  fdled  honeycomb  tend  to  be  complicated  and  are  valid  only 
over  a  limited  range  of  infdl  moduli,  often  researchers  default  to  finite  element  modeling  when 
investigating  such  structures.  Such  methods,  however,  have  limited  use  in  optimization.  Works 
involving  the  modeling  and  experimental  and  FEA  validation  of  filled  and  unfilled  honeycomb 
and  other  cellular  structures  include  those  by:  [42,  43,  30,  44,  45,  46,  47,  48,  49,  50]  [51] 

Regarding  the  modeling  of  filled  honeycombs,  more  recent  works  by  Grarniiller,  Pagitz,  and  Vos 
and  their  colleagues  have  resulted  in  filled  and  pressurized  honeycomb  models.  In  particular, 
Pagitz’  work,  inspired  by  the  nastic  movements  of  plants,  utilizes  pressure  differences  between 
honeycomb  cells  to  achieve  deformation.  Their  model  assumes  honeycomb  made  of  rigid  beams 
connected  with  springs.  The  properties  of  the  structure  and  deformation  are  then  calculated  for 
various  internal  pressures.  [52]  Grarniiller  et  al.  recently  used  this  method  to  demonstrate  a 
morphing  airfoil  concept.  [53]  Similar  models  developed  by  Vos  and  Barrett  relating  the  internal 
pressure  of  a  honeycomb  cell  to  mechanical  work  have  also  been  used  in  the  topology 
optimization  of  actuators  for  morphing  leading  and  trailing  edge  flaps.  Similar  to  Pagitz’  work, 
changes  in  pressure  in  the  honeycombs  result  in  large  scale  deformation  of  the  wing.  Vos  and 
Barrett  showed  that  this  style  of  analytical  model  correlates  well  with  experimental  results  and 
FEA  analysis.  [54,  55,  56,  57]  The  major  difference  of  these  methods  with  respect  to  the 
proposed  application  is  the  forces  exerted  on  the  honeycomb  walls  by  internal  pressure 
differences  are  uniform  and  normal  to  the  cell  walls,  neither  of  which  hold  true  for  solid  infdl. 

To  aid  in  the  design  and  optimization  of  a  cellular  skin,  it  is  helpful  to  develop  an  accurate 
analytical  model  of  the  system.  The  model  can  then  be  used  to  calculate  in-plane  composite 
mechanical  properties  such  as  Young’s  modulus  and  Poisson’s  ratio  which  in  turn  can  be  used  by 
optimization  routines.  Understanding  the  mechanics  of  honeycomb  and  the  methods  used  to 
model  them  are  essential  in  developing  an  improved  model  of  fdled  honeycomb.  For  this  reason 
a  brief  description  of  three  empty  honeycomb  models  foUows.  The  three  models  discussed  here 
predicting  the  mechanical  properties  of  unfdled  honeycomb  are  listed  in  Table  10  (the 
coefficients  of  which  are  to  be  inserted  into  Equation  (27)  with  B2=D2=0)  with  variable 
definitions  shown  in  Figure  5 1 .  The  first,  by  Masters  and  Evans,  calculates  the  deflection  of  a 
honeycomb  ceU  by  first  calculating  deflection,  axial  strain,  and  shear  using  beam  theory 
separately  for  each  member  before  combining  them  into  a  single  equation.  [58]  This  results  in 
the  lowest  predicted  Young’s  modulus  (see  Table  7). 

The  second  model,  the  simplified  version  by  Gibson  and  Ashby,  is  the  most  widely  used  due  to 
its  relative  accuracy  and  simplicity.  [59]  It  is  similar  to  more  complicated  models,  except  that  it 
neglects  all  deformation  besides  bending.  Although  simplified,  it  is  relatively  accurate  and 
consistent  with  other  models.  The  third  model  presented  in  Table  10  is  the  unabridged  version  by 
Gibson  and  Ashby  allowing  for  shear  deformation,  which  has  been  successfully  used  as  a  basis 
by  Olympic  and  Gandhi  to  model  hybrid  and  accordion  cellular  cores.  [25,  26,  27] 
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While  there  are  many  models  that  predict  the  properties  of  unfdled  honeycomb,  see  Table  10, 
fewer  predict  the  properties  of  filled  honeycomb.  Moreover,  those  that  do  are  only  accurate  at 
very  low  infdl  moduli.  Equation  (27)  is  a  generic  model  used  to  calculate  the  effective  stiffness 
of  fdled  honeycomb  composite  by  separating  the  contributions  of  the  honeycomb  and  the  infdl. 


^cy  ~  A  A 


(27) 


In  Equation  (27)  E  is  Young’s  modulus,  B  and  D  are  geometric  coefficients  listed  in  Table  10, 
subscript  c  represents  composite  quantities,  h  represents  honeycomb  quantities,  i  represents  infill 
quantities,  and  x  and  y  are  coordinate  directions.  Note  that  this  equation  represents  a  generic 
model  that  captures  the  separate  contribution  of  the  infill  and  the  honeycomb  to  the  modulus  of 
the  combined  system.  This  is  valid  for  rule-of-mixtures  models  as  well  as  more  complicated 
models  as  described  below. 


Of  particular  interest  to  this  work  is  the  honeycomb  mechanics  model  developed  by  Abd  El- 
Sayed  for  honeycombs  with  a  solid  infill.  [60,  61]  Their  work  is  the  earliest  presented  and  the 
bases  for  many  other  honeycomb  models,  including  the  often  cited  models  by  Gibson  and  Ashby. 
The  model  uses  Castigliano’s  method  to  calculate  the  internal  strain  energy  associated  with  each 
honeycomb  member  similarly  to  Gibson  and  Ashby’s  simplified  model  however  includes  axial 
deformation.  Coefficients  accounting  for  the  contributions  of  the  honeycomb  and  infill  when 
modelling  filled  honeycomb  are  listed  in  Equation  (28).  El-Sayed’s  model  predicts  Young’s 
modulus  in  two  directions  as  well  as  Poisson’s  ratio  in  the  elastic,  elasto-plastic,  and  plastic 
deformation  regions.  The  model  is  shown  to  correlate  well  with  experimental  data  with  an  infill 
modulus  between  zero  and  four  orders  of  magnitude  or  more  lower  than  that  of  the  honeycomb. 
The  equations  shown  in  Equations  (28)  through  (31)  were  re-derived  using  the  methods  of  Abd 
El-Sayed  due  to  typos  in  the  original  articles. 
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El-Sayed’s  model  provides  acceptable  approximations  of  Ecx  and  Ecy  when  Ei  is  much  less  than 
Eh-  When  Ei  is  relatively  large,  Equation  (27)  greatly  overestimates  the  contribution  of  the  infill 
resulting  in  overestimates  of  Ecx  and  Ecy.  This  is  due  primarily  to  the  forces  exerted  on  the  beam 
members  by  the  infill  being  non-uniformly  distributed  along  the  lengths  of  the  honeycomb 
members  and  neglecting  the  resulting  mechanical  advantage.  For  this  reason,  a  new  model  was 
developed  in  the  same  manner  as  the  El-Sayed  model  to  more  accurately  represent  the  internal 
distributed  forces  resulting  from  deformation  of  a  fdled  honeycomb  structure  and  improve  the 
estimate  of  composite  modulus  for  large  values  of  infdl  modulus. 

Table  10:  Current  Unfilled  Honeycomb  Model  Coefficients 
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3.6.2  Anal5d;ic  Filled  Honeycomb  Model  Derivation 

If  the  honeycomb  is  viewed  as  a  system  comprised  of  beam  elements  rigidly  connected  at  each 
end,  significant  mechanical  advantage  is  generated  by  the  infill  resisting  deformation  of  the 
honeycomb  walls  via  a  distributed  load.  Forces  exerted  by  the  infill  due  to  strain  near 
connection  points  have  a  much  lower  impact  on  the  effective  stiffness  of  the  composite  than 
forces  far  from  connection  points.  It  is  postulated  that  neglecting  this  mechanical  advantage  is 
partially  the  cause  of  the  limitation  of  low  infill  moduli  of  available  models.  To  account  for  this 
mechanical  advantage,  the  presented  model  estimates  the  effect  of  the  infill  by  applying 
distributed  loads  directly  to  honeycomb  beam  members  as  functions  of  strain.  Calculating  strain 
energy  in  the  infill  directly  is  then  unnecessary. 

The  following  model  predicting  the  effective  in-plane  Young’s  moduli  of  filled  honeycomb  is 
derived  using  thin  beam  theory  and  energy  methods.  Materials  are  assumed  to  be  linearly  elastic 
with  each  honeycomb  member  within  the  unit  cell,  see  Figure  51,  considered  to  be  a  slender 
beam  rigidly  attached  at  the  ends. 
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Figure  51:  Honeycomb  and  Unit  Cell  Geometry  and  Dimensions 

The  majority  of  the  external  load  is  assumed  to  be  carried  by  the  honeycomb  and  Castigliano’s 
method  is  used  to  calculate  the  force-deflection  relationship  of  the  honeycomb  in  each  desired 
direction: 
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where  S  is  the  deflection  in  the  direction  and  at  the  point  of  the  applied  force/  m  is  the  number 
of  members,  is  the  axial  load,  E  is  Young’s  modulus,  A  is  the  cross-sectional  area  of  the 
member,  M  is  the  moment,  /  is  the  second  moment  of  area,  A:  is  a  correction  coefficient,  V  is  the 
shear  force,  G  is  the  shear  modulus,  T  is  the  torque,  and  J  is  the  polar  moment  of  inertia.  The 
subscript  j  represents  the  direction  of  the  deflection  being  sought;  in  this  case  either  0,  V2,  0,  or 
0+V2  shown  as  a  in  Figure  51. 


12  (33) 

A  =  cd 

With  the  unit  cell  corresponding  to  an  a  of  0  shown  in  Figure  5 1,  the  forces  on  the  slanted  and 
axial  honeycomb  segments  within  the  unit  cell  are  similar  to  those  shown  in  Figure  52  and 
Figure  53. 
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Figure  52:  Forces  on  Slanted  Members 


y,r 

-h 


Figure  53:  Forces  on  Horizontal  Members 


The  subscript  i  indicates  a  force  generated  by  straining  the  infill  and  the  subscript  0  refers  to  an 
applied  force.  The  moment  mg  is  found  such  that  a  fixed  end  condition  is  imposed.  Although 
some  SMPs  can  exhibit  nonlinear  behavior,  for  simplicity  the  infill  material  is  assumed  linear 
here  thus  limiting  the  application  of  the  following  equations  to  small  deflections  within  the  linear 
elastic  region  of  both  the  honeycomb  and  infill  materials.  At  large  strains,  the  presented 
equations  become  increasingly  erroneous  due  to  nonlinear  effects  and  changing  honeycomb 
geometry.  The  forces  due  to  the  infill  in  the  a=0  direction  are 
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The  subscript  i  indicates  a  force  due  to  the  infill,  x  (a  =  0)  or  y  (a  =  n/2)  indicates  the  direction 
of  the  force,  a  or  I  indicates  the  type  of  member  the  force  acts  on,  either  slanted  or  horizontal 
(see  Figure  51).  The  remaining  quantities  are 
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It  should  be  noted  that  the  infdl  also  exerts  shear  and  torque  forces  locally  on  the  honeycomb 
members  as  well  as  the  axial  forces  listed  in  Equation  (34).  These  forces  were  found  to  have 
little  influence  on  the  effective  stiffness  of  the  composite  and  are  thus  neglected  to  simplify  the 
equations.  Decomposing  the  forces  in  Equation  (34)  into  components  in  the  r  and  5  directions 
and  neglecting  torque  and  shear  from  Equation  (32),  the  resulting  deflection  equation  in  the  a=0 
direction  is 


^0  if 0,0s  ) ^^/Ir  (fo,Or  ’  fo,iyar  ’  f  0,vcr  ’  fo,iynr  )  if  0,0s  ’  fo,iyas  ’  fo,iyl\s  ’  fo,ixs  ) 


(36) 


The  subscripts  a,  1, 11,  and  12  indicate  the  respective  honeycomb  members  within  the  unit  cell 
while  the  subscripts  r  and  s  indicate  direction  as  shown  in  Figure  52  and  Figure  53.  In  Equation 
(36)  symmetry  yields  5/;  and  5/2  are  equivalent;  thus,  only  one  needs  to  be  solved  for  and  a  factor 
of  2  applied.  The  lower  case  forces,/,  correspond  to  the  decomposed  forces  in  Equations  (34)  in 
the  r  and  s  coordinate  system.  For  example,  the  force  deflection  equation  of  a  slanted  beam  due 
to  the  infill  in  the  x  direction  (Fojx)  is  decomposed  based  on  the  honeycomb  geometry  (0)  into 
fo,ixr  and  fo_ixs-  Inserting  the  decomposed  forces  into  Equation  (36)  and  the  rearranging  terms 
yields  the  total  deflection  of  the  unit  cell  in  the  a=0  direction: 
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Notice  that  Equation  (37)  is  not  solved  for  the  deflection,  5o,  explicitly  for  clarity  since  the  forces 
due  to  the  infill  are  dependent  upon  deflection.  Using  Hooke’s  law.  Equation  (37),  and  the 
dimensions  of  the  unit  cell  the  tensile  modulus  in  the  a=0  direction  is 
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If  the  effect  of  the  infdl  is  ignored.  Equation  (38)  can  be  simplified  and  shown  to  be  equivalent 
to  that  of  Abd  El-Sayed  shown  in  Table  10  and  Equation  (27).  This  method  is  then  used  to  solve 
for  the  effective  stiffness  of  the  filled  honeycomb  in  the  remaining  three  directions:  a  =  {%,  0, 
and  0+V2},  the  equations  for  which  are  shown  in  the  appendix.  A  spline  interpolation  between 
each  modulus  is  then  used  to  estimate  the  modulus  of  the  honeycomb/SMP  composite  in  any 
direction.  If  the  infill  modulus  is  set  to  zero,  the  effective  in-plane  stiffness  of  empty  honeycomb 
is  captured. 
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A  similar  process  can  be  employed  to  solve  for  the  effeetive  shear  modulus  of  the  eomposite. 
The  foree  due  to  the  infill  is 


To 


The  foree  deflection  equation  for  shear  in  the  a  =  (0,%)  direction  is 
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Decomposing  Equation  (39)  into  r  and  s  components,  inserting  them  into  Equation  (40),  and 
rearranging  the  terms  yields: 
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(41) 


Using  the  definition  of  shear  modulus,  Equation  (41),  and  the  dimensions  of  the  unit  cell  the 
effeetive  shear  modulus  is  then 


G 


-^0,0  2^0 
Aq  c{a  +  Xq  ) 


(42) 


This  method  can  also  be  used  to  calculate  the  effeetive  shear  modulus  of  the  filled  honeycomb 
composite  in  the  remaining  three  directions  a  =  {(V2,0),  (0,0+%),  and  (0,0+%)},  shown  in  the 
appendix.  A  spline  interpolation  between  eaeh  point  is  then  utilized  to  predict  the  shear  modulus 
at  any  angle.  If  the  infill  modulus  is  set  to  zero,  the  effeetive  in-plane  shear  modulus  of  empty 
honeyeomb  is  found. 


3.6.3  Analytic  Filled  Honeycomb  Model  Results 

Table  1 1  and  Figure  54  show  the  results  of  the  presented  analytical  model  predicting  Young’s 
moduli  of  empty  honeyeomb  eompared  to  the  experimentally  determined  values  at  angles  a  =  0, 
0,  0+"/2,  and  %.  Analytieal  results  are  indieated  by  blue  dots,  a  spline  interpolation  between  data 
points  is  represented  by  a  dotted  blue  line,  and  experimental  values  and  their  associated  error 
bars  are  shown  in  red  in  Figure  54.  For  Figure  54  through  Figure  58  the  radial  axis  is  Young’s  or 
shear  modulus  in  MPa  and  the  polar  axis  is  a  from  Figure  5 1  in  radians. 
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Figure  54:  Empty  Honeycomb  Young’s  Modulus  (MPa) 
Table  11:  Empty  Honeycomb  Young’s  Modulus 


Direction 

(a) 

Experimental 

(kPa) 

Analytical 

(kPa) 

0 

55.5±11.7 

26.4 

71/ 

/l2 

472±58 

206 

Til 

/6 

783±130 

511 

0 

637±72 

658 

71/ 

12 

20.1±1.8 

8.95 

0+^2 

NA 

667 

0  =  0.799  radians 


Both  the  model  and  experimental  data  show  a  signifieant  increase  in  stiffness  in  directions 
approaching  the  0  direction  of  the  honeycomb.  At  an  angle  of  a  =  0,  the  applied  load  is  axial  to 
one  of  the  slanted  members  of  the  honeycomb  unit  cell.  Since  the  slanted  members,  /,  are 
significantly  longer  than  the  horizontal  members,  a,  and  the  dominant  mode  of  deformation  is 
bending  of  these  members,  total  deformation  is  at  a  minimum  when  one  of  the  slanted  members 
is  loaded  axially,  resulting  in  increased  stiffness.  If  a  standard  honeycomb  where  the  length  of 
member  a  equals  the  length  of  members  I  and  a  =  60°  is  modeled,  the  predicted  stiffness  is 
nearly  isotropic  in-plane.  From  Figure  54  the  maximum  stiffness  of  the  empty  honeycomb  was 
found  to  be  when  a  =  30°  experimentally  whereas  the  model  predicts  the  maximum  modulus 
should  occur  at  a  =  0.  There  are  two  possible  contributing  factors  for  this  discrepancy.  The 
spline  interpolation  between  data  points  (the  dotted  blue  line  in  Figure  54)  may  be  inappropriate 
and  misleading  resulting  in  a  more  pointed  curve.  Also,  due  to  the  geometry  of  the  honeycomb 
there  are  large  errors  when  measuring  the  modulus  of  the  material  at  arbitrary  angles  as  indicated 
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by  the  large  error  bars  in  Figure  54.  It  is  possible  that  this  discrepancy  also  occurs  when  the 
infill  modulus  is  nonzero;  however  no  experimental  data  is  available. 


Equations  (39)  through  (42)  are  then  used  to  calculate  the  shear  modulus  of  empty  honeycomb 
by  setting  the  infill  modulus  to  zero.  From  Table  12,  the  shear  modulus  of  honeycomb  follows  a 
similar  directional  pattern  as  Young’s  modulus,  however  less  pronounced. 

Table  12:  Empty  Honeycomb  Shear  Modulus 

Direction  Analytical 
(g) _ (kPa) 

0,"/2  8.48 

0,0+%  7.63 

%,0  25.4 

0+%,0  7.58 


Table  13  and  Figure  55  show  experimental  and  analytical  Young’s  moduli  predictions  for  filled 
honeycomb  above  Tg.  The  range  of  values  of  stiffness  in  various  directions  is  significantly 
reduced  compared  to  empty  honeycomb  as  the  addition  of  SMP  has  an  averaging  effect  on  the 
stiffness  of  the  composite.  The  difference  between  the  analytically  predicted  Young’s  modulus 
in  the  a=0  and  %  directions  and  the  experimentally  determined  values  are  less  than  6.5  and  4.3 
percent,  respectively. 


Table  13:  Composite  Young’s  Modulus  >Tg 
Direction  (a)  Experimental  (MPa)  Analytical  (MPa)  El-Sayed  (MPa) 


0 

33.9±4.1 

36.1 

72.74 

0 

NA 

20.9 

NA 

% 

11.8+4.2 

12.3 

24.80 

0+% 

NA 

21.2 

NA 

0.799  radians  for  the  presented  case 
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Figure  55:  Composite  Young’s  Modulus  >Tg  (MPa) 

While  the  models  available  in  the  literature  ealculating  Young’s  modulus  of  empty  honeyeomb 
are  very  similar  to  the  presented  model,  the  differenees  become  apparent  when  the  infill  is 
considered.  In  Table  13,  the  filled  honeycomb  model  by  El-Sayed  predicts  more  than  double  the 
experimentally  determined  composite  stiffness  in  both  the  a  =  0  and  nil  directions. 

When  the  shear  modulus  of  the  composite  above  Tg  is  graphed  in  a  similar  manner,  the  familiar 
flower  shape  returns.  The  error  between  the  experimentally  determined  moduli  and  the 
presented  model  are  2.9%  in  the  (0,"/2)  direction  and  12.3%  in  the  ("/2,0)  direction.  The  order  of 
magnitude  difference  in  the  shear  modulus  of  the  material  in  the  (0,%)  and  ("/2,0)  directions 
compared  to  the  (0,0+"/2)  and  (0+V2,0)  directions  is  believed  to  be  due  to  the  construction  of  the 
honeycomb.  Honeycomb  members  a  are  twice  as  thick  as  members  /,  greatly  increasing  bending 
stiffness  in  certain  directions. 


Table  14:  Composite  Shear  Modulus  >Tg 


Direction  (a) 

Experimental  (MPa) 

Analytical  (MPa) 

0,% 

13.9±3.5 

14.3 

0,0+% 

NA 

0.230 

%,o 

13.0+2.7 

11.4 

0+%,0 

NA 

4.75 
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Figure  56:  Composite  Shear  Modulus  >Tg  (MPa) 

In  the  cold  state,  the  stiffness  of  the  polymer  is  within  1.5  orders  of  magnitude  of  the  aluminum 
honeycomb.  When  the  two  materials  are  this  similar  in  moduli  the  dominant  mode  of 
deformation  is  no  longer  bending  and  the  presented  model  becomes  invalid.  This  is  clearly 
shown  in  Table  15  and  Figure  57,  where  the  circular  red  dashed  line  represents  the  stiffness  of 
the  composite  calculated  using  rule  of  mixtures.  As  the  stiffness  of  the  polymer  approaches  the 
stiffness  of  the  honeycomb  the  composite’s  stiffness  approaches  that  of  a  solid  isotropic  flat  plate 
and  rule  of  mixtures  becomes  more  accurate  than  the  presented  model.  In  the  presented  case 
below  Tg  rule  of  mixtures  yields  1.86E9  Pa. 

Table  15:  Composite  Young’s  Modulus  <Tg 


Direction  (a) 

Experimental  (MPa) 

Analytical  (MPa) 

El-Sayed  (MPa) 

0 

2190±110 

2480 

4987 

0 

NA 

1390 

NA 

% 

2040±140 

846 

1700 

e+% 

NA 

1420 

NA 
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Figure  57:  Composite  Young’s  Modulus  <Tg  (MPa) 

As  with  Young’s  modulus,  the  error  between  the  presented  model  and  experimentally 
determined  shear  moduli  values  is  much  greater  when  the  composite  is  below  Tg,  see  Figure  58. 
Without  further  experimental  investigations  of  the  stiffness  of  the  composite  at  angles  other  than 
0  and  however,  it  is  uncertain  if  the  predicted  shapes  in  Figure  54  tWugh  Figure  58  are 
realistic. 
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3.6.4  Honeycomb  Geometry  Optimization 


Using  El-Sayed's  equations  for  the  equivalent  Young’s  moduli  and  flat  plate  theory,  the  out  of 
plane  deflection  of  the  honeycomb  panel  can  be  derived.  The  weight  of  the  system  is  derived  by 
computing  the  total  number  of  cells  in  a  panel,  N. 


N  = 


A 

2yoia  +  Xo) 


(43) 


where  A  is  the  area  of  the  entire  panel.  The  area  of  the  infill  and  honeycomb  in  a  unite  cell  is 
then: 


=  2(yo  -  d)(Xo  -  d)  +  2((a  -  d)(yo  -  d))  (44) 

Ah  =  (a  +  Xo)  -  Ai  (45) 

The  density  of  the  infill  and  honeycomb  are  then  used  to  calculate  the  weight  (IT)  of  the  skin. 

W  =  (A,p,+AhPh)cN  (46) 


3.6.4. 1  Optimization 

The  honeycomb  unit  cell  geometry  was  optimized  to  minimize  the  Young’s  modulus  in  the  x  and 
y  direction.  This  will  ultimately  reduce  the  actuation  energy  needed  to  morph  the  composite 
skin.  The  panel  of  interest  was  based  on  the  MAS  Phase  II  design  requirements  with  dimensions 
of  .508  X  .381  m,  constrained  to  deflect  2.5  mm  in  the  middle  of  the  plate  when  a  19,152  N/m^ 
load  is  applied.  The  infdl  SMP  has  a  modulus  of  1.6e9  Pa,  the  soft  state  of  the  polymer,  while 
the  honeycomb  has  a  modulus  of  73.1  GPa. 

The  design  variables  are  the  horizontal  wall  length,  a,  slanted  wall  length,  /,  wall  thickness  d,  and 
cell  angle,  0,  see  Figure  51.  The  depth  of  the  cell  and  skin  weight  were  not  included  as  design 
variables  as  only  in  plane  performance  was  considered. 

The  objective  of  the  design  is  to  minimize  Equation  (47)  resulting  in  the  lowest  actuation  cost  in 
both  the  X  andy  direction. 


Ec  H"  E^ 

Cx  Cy 

(47) 

The  cell  angle  was  constrained  to  en  sure  the  optimized  design  would  be  capable  of  meeting  the 
MAS  strain  requirements.  The  cell  wall  thickness,  d,  had  a  lower  bound  constraint  of  0.2  mm 

which  is  the  lowest  manufactured  honeycomb  cell  wall  width, 
lengths  were  also  given  an  arbitrary  upper  bound. 

The  horizontal  and  slanted  wall 

55°  <  e»  <  75° 

(48) 

0.2  mm  <  d 

(49) 

l,a  <  100  mm 

(50) 
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The  weight  of  the  panel  was  initially  constrained  by  using  the  depth  of  the  cell  as  a  design 
variable,  however  it  was  discovered  that  both  the  weight  and  deflection  constraints  could  not  be 
simultaneously  met  with  the  infill  in  the  soft  state.  The  design  with  the  lowest  weight 
maintaining  a  panel  center  point  deflection  below  2.5  mm  weighed  8.3  kg;  significantly  violating 
the  MAS  requirement. 

The  optimization  was  run  with  MatLAB©  Optimtool.  Optimtool  selects  the  ‘best’  algorithm  to 
use  which  was  a  mixture  of  the  sequential  quadratic  programming,  Quasi-Newton,  and  line 
search  algorithms. 


3. 6.4.2  Optimization  Results 

A  series  of  sensitivity  studies  were  conducted  by  plotting  Exc  and  Eye  verses  a  set  of  two  design 
variables.  Figure  59  and  Figure  60  look  at  the  effect  wall  thickness,  d,  and  slanted  wall  length,  /, 
of  the  honeycomb  have  on  the  stiffness  of  the  honeycomb  in  the  x  and  y  directions,  Exc  and  Eye. 
As  expected  since  the  major  mode  of  deformation  within  the  honeycomb  is  bending,  the  wall 
thickness,  d,  has  a  large  impact  on  the  stiffness  of  the  composite  in  both  directions.  Still 
significant  however  less  drastic,  shorter  slanted  honeycomb  members  result  in  higher  stiffnesses 
in  the  x  direction  and  at  high  values  of  d  in  the  y  direction. 


I  -  Slanted  Wall  Length,  mm 


d  -  Wall  Thickness,  mm 


Figure  59:  Young’s  modulus  in  the  x  direction,  Exc,  of  the  optimal  cell  as  a  function  of  the 

slanted  wall  length,  I,  and  cell  wall  thickness,  d 


I  -  Slanted  Wall  Length,  mm 


d  -  Wall  Thickness,  mm 


Figure  60:  Young’s  modulus  in  the_v  direction.  Eye,  of  the  optimal  cell  as  a  function  of  the 
slanted  wall  length,  /,  and  the  cell  wall  thickness,  d 
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Figure  61  and  Figure  62  show  the  effect  of  the  horizontal  wall  length,  a,  with  respect  to  wall 
thickness,  d,  on  the  stiffness  of  the  composite  in  the  x  and  y  directions  respectively.  Opposite  to 
the  effect  of  the  slanted  wall  length,  shorter  horizontal  members  result  in  less  stiff  geometries  in 
the  X  direction.  In  the  y  direction,  shorter  horizontal  members  result  in  a  stiffer  composite. 


10^ 


a  -  Horizontal  Wall  Length,  mm 


d  -  Wall  Thickness,  mm 


Figure  61:  Young's  modulus  in  the  x  direction,  Exc-,  of  the  optimal  cell  as  a  function  of  the 
horizontal  wall  length,  a,  and  cell  wall  thickness,  d 


d  -  Wall  Thickness,  mm 

a  -  Horizontal  Wall  Length,  mm 


Figure  62:  Young’s  modulus  in  the  j  direction.  Eye,  of  the  optimal  cell  as  a  function  of 
horizontal  wall  length,  a,  and  cell  wall  thickness,  d 

Figure  63  and  Figure  64  show  the  importance  of  honeycomb  lengths  a  and  I  on  Exc  and  Eye.  In 
both  directions,  the  length  of  the  slanted  member,  /,  is  significantly  more  influential  over  the 
effective  stiffness  of  the  honeycomb  than  the  horizontal  member,  a. 
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Youngs  Modulus  in  x  Direction 


Figure  63:  Young's  modulus  in  the  jc  direction,  Exc,  of  the  optimal  cell  as  a  function  of  the 

slanted  wall  length,  I,  and  cell  angle,  6 


Youngs  Modulus  in  y  Direction 


I  -  Slanted  Wall  Length,  mm 


a  -  Horizontal  Wall  Length,  mm 


Figure  64:  Young's  modulus  in  the  j  direction,  Exc,  of  the  optimal  cell  as  a  function  of  the 

slanted  wall  length,  I,  and  cell  angle,  G 

Over  the  range  of  angles  investigated  (55°<0<75°),  the  internal  honeycomb  angle,  0,  has  little 
effect  on  the  stiffness  in  the  direction,  Exc,  shown  in  Figure  65  and  Figure  67.  In  the  y  direction, 
however,  the  angle  is  much  more  significant  with  larger  angles  resulting  in  higher  stiffness  as 
shown  in  Figure  66  and  Figure  68. 


Figure  65:  Young’s  modulus  in  the  x  direction,  Exc,  of  the  optimal  cell  as  a  function  of 

horizontal  wall  length,  a,  and  cell  angle, 
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a  -  Horizontal  Wall  Length,  mm 


0  -  Cell  Angle,  deg 


Figure  66:  Young’s  modulus  in  the  j  direction,  Eye,  of  the  optimal  as  a  function  of  the 

horizontal  wall  length,  a,  and  cell  angle, 


I  -  Slanted  Wall  Length,  mm 


0  -  Cell  Angle,  deg 


Figure  67:  Young’s  modulus  in  the  x  direction,  Exc,  with  the  slanted  wall  length,  /,  and  cell 

angle, 


Figure  68:  Young's  modulus  in  the_v  direction.  Eye,  with  the  slanted  wall  length,  /,  and  cell 

angle, 
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80 


10 


0  -  Cell  Angle 


d  -  Wall  Thickness,  mm 


Figure  69:  Young’s  modulus  in  the  x  direction,  Exc,  of  the  optimal  cell  as  a  function  of  the 

cell  wall  thickness,  d,  and  cell  angle, 


0  -  Cell  Angle 


d  -  Wall  Thickness,  mm 


Figure  70:  Young’s  modulus  in  the  j  direction.  Eye,  of  the  optimal  cell  as  a  function  of  the 

cell  wall  thickness,  d,  and  cell  angle, 

The  sensitivity  study  illustrates  that  the  largest  contributing  factor  to  the  overall  effective  moduli 
of  the  composite  is  the  wall  thickness.  Large  angles  of  <5^  caused  Eycto  increase  and  Exc  to 
decrease  as  expected. 

The  optimization  was  started  at  various  initial  starting  points  throughout  the  range  of  each 
variable  producing  different  /,  d,  and  a  values  depending  on  initial  conditions.  The  angle  ^^'was 
consistently  limited  by  the  upper  bound  constraint.  The  ratio  of  a:l  and  a:d  were  also  constant 
for  each  solution  at  0.3  and  26.4  respectively.  Presented  in  Table  16,  are  two  sets  of  results. 
Note  that  the  moduli  and  d  are  the  same. 


Table  16:  Optimal  Design  Points 


Exc  Eye  ^  I  d  a 

GPa  GPa  deg  mm  mm  mm 


1.55  1.36  75  16.19  2.01  52.14 

1.55  1.36  75  1.619  .201  5.28 
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3.7  System  Design  and  Optimization 

3.7.1  Heating  Patterns 

Honeycomb  stractures  are  inherently  flexible  both  in-plane  and  in  bending  due  to  their  generally 
thin  walled  construction  when  not  paired  with  face  sheets.  The  structure  can  be  made  relatively 
stiff,  however,  by  filling  the  voids  in  the  honeycomb  with  a  rigid  material  effectively  resulting  in 
a  solid  sheet.  If  the  honeycomb  is  filled  with  a  thermally  activated  SMP  material  and  each  cell 
heated  individually,  thus  providing  cellular  level  stiffness  control,  the  filled  honeycomb 
composite  becomes  a  stiffness  controllable,  reconfigurable  skin  system.  Patterns  of  heated  cells, 
such  as  those  shown  in  Figure  71,  can  then  be  used  to  alter  the  macroscopic  properties  of  the 
material. 


Pattern  5 


Pattern  6 

Figure  71.  Honeycomb  Heating  Patterns 


Pattern  7 


3.7.2  FEA  Modeling 

To  estimate  the  effective  stiffness  of  various  heating  patterns  and  assess  their  capabilities,  a 
model  that  is  eight  cells  wide  by  five  cells  high  was  generated  in  Abaqus  using  the  geometry  of 
the  experimentally  characterized  honeycomb  detailed  previously.  The  aluminum  honeycomb 
was  modelled  using  beam  elements  while  the  SMP  infill  was  modelled  using  2D  shell  elements. 
All  proper  convergence  studies  were  performed  ensuring  the  model  size,  boundary  conditions, 
and  mesh  densities  were  adequate.  Each  of  the  seven  heating  patterns  shown  in  Figure  71  was 
modelled  as  well  as  empty  honeycomb.  Figure  72  is  representative  of  the  data  obtained  showing 
Pattern  7. 
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I  Step:  Step-1,  Stepl 

Ificrement  liStepTlme=  1.000 
Primery  Van  S,  Mises 

Deformed  Van  U  Deformation  Scale  Factor:  +4.704e  +  00 


Figure  72:  FEA  Mises  Stress  Results  for  Pattern  7  (Figure  71)  Under  Shear 

Numerical  values  for  the  effective  stiffness  in  the  x  and  y  directions  for  epoxy  and  PEEK  filled 
aluminum  honeycomb  are  listed  in  Table  17  and  shown  in  Figure  73.  Material  properties  for  the 
two  polymers  can  be  found  in  Table  18.  Notably,  for  heating  pattern  three,  epoxy  produces  a 
composite  stiffer  in  the  x  direction  than  PEEK  but  less  stiff  in  the  y  direction  with  the  reverse 
true  for  Pattern  4.  For  Pattern  6,  epoxy  and  PEEK  result  in  nearly  identical  moduli  in  the  x 
direction.  In  the  y  direction  however,  PEEK  produces  a  modulus  less  than  half  that  of  epoxy. 
Comparing  the  stiffnesses  produced  by  these  two  polymers  clearly  illustrates  the  complexity  of 
material  selection  for  the  proposed  system.  The  benefit  of  variable  heating  patterns  is  also 
clearly  demonstrated  when  comparing  the  stiffnesses  of  patterns  3  and  4.  As  expected.  Pattern  3 
is  stiff  in  the  y  direction  and  soft  in  the  x  while  Pattern  4  is  the  opposite. 

Table  17:  FEA  Predicted  Effective  Young’s  Modulus  for  PEEK  and  Epoxy  Infill 

_ 1 _ 2 _ 3 _ 4 _ 5 _ 6 _ 7 _ 8 

Ex  Epoxy  1.32E-F09  2.53E-F07  9.92E-F07  5.00E-F08  1.05E-F08  1.01E-F08  1.57E-F08  1.55E4-05 

_ PEEK  3.61E-F09  6.05E-F06  6.69E-F07  1.31E-F09  1.21E-F08  1.05E-F08  3.09E-F08  1.55E4-05 

Ey  Epoxy  1.31E-F09  2.01E-F07  4.46E-F08  7.98E-F07  4.87E-F07  4.93E-F07  6.63E-F07  5.30E4-04 

PEEK  3.60E-F09  2.94E-F06  1.17E-F09  3.72E-F07  2.61E-F07  2.38E-F07  8.50E-F07  5.30E4-04 
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Young's  Modulus  (Pa) 


Young's  Modulus  (Pa) 


X  Direction  X  Direction 


Figure  73:  FEA  Predicted  Young’s  Modulus  for  Each  Heating  Pattern  in  the  X  and  Y 
Directions  for  Epoxy  (Left)  and  PEEK  (Right)  Eilled  Aluminum  Honeycomb 

3.7.3  Material  Selection 


With  weight  being  the  most  limiting  factor  in  the  design,  material  specific  strength  is  an 
important  material  property.  Several  shape  memory  polymers  as  well  as  silicone  rubber  are 
listed  in  Table  18  with  their  associated  transition  temperatures,  densities,  stiffnesses,  and  specific 
strengths.  The  highest  listed  specific  strength  in  the  cold  state  is  PEEK  indicating  this  material 
could  result  in  the  lowest  weight.  CRG’s  Veriflex  and  VeriflexE  have  the  lowest  specific 
strength  above  Tg  indicating  their  potential  to  have  the  lowest  actuator  power  requirements.  As 
indicated  by  the  FEA  analysis  and  skin  optimization,  material  selection  is  complexly  related  to 
desired  heating  patterns,  which  are  determined  by  desired  performance  characteristics  of  the 
skin.  The  only  SMP  listed  that  meets  Altera  Corporation’s  temperature  rating  for  military  grade 
devices,  Table  19,  is  PEEK.  Other  factors  not  included  in  this  study  that  would  need  to  be 
considered  are  cost,  manufacturability,  reliability,  maintenance,  toxicity,  failure  modes,  and 
environmental  factors  such  as  resistance  to  dirt,  oxidation,  and  moisture. 

Table  18:  Specific  strength  of  available  SMPs 


<Tg ^ 


Name 

T 

°c 

P 

kg/m^ 

Modulus 

Pa 

Specific 
Strength 
MPa  m^/kg 

Modulus 

Pa 

Specific 
Strength 
MPa  m^/kg 

PEEK  polyether  ether  ketone 

143 

1320 

3.60E9 

2.73 

2.60E6 

0.00197 

CRG  VeriflexE  epoxy 

100 

1153 

3.03E9 

2.63 

6.60E5 

0.00057 

MHI MM4510  urethane-ester 

45 

1200 

3.00E9 

2.50 

4.50E6 

0.00375 

EPON  826/Jeffamine  D230 

90 

1153 

2.70E9 

2.34 

2.00E7 

0.01735 

polyethylene  terephthalate  (PET) 

67-81 

1380 

2.8E9-4.5E9 

2.029-3.261 

5.00E8 

0.36232 

MHI  MM4520  urethane-ether 

45 

1250 

2.00E9 

1.60 

4.00E6 

0.00320 

polyurethane  MM3520 

36.7 

1250 

1.60E9 

1.28 

2.30E6 

0.00184 

CRG  Veriflex  polystyrene 

62 

920 

1.14E9 

1.24 

4.40E5 

0.00048 

EPON  826/Jeffamine  D230 

78 

1153 

1.31E9 

1.14 

1.90E7 

0.01651 

silicone  rubber 

NA 

1100-2300 

1.0E6-5.0E7  0.00044-0.045 

1.0E6-5.0E7 

0.00043-0.045 
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Table  19:  Altera  Corporation’s  temperature  requirements  for  military  grade 

programmable  logic  devices 


Industry  Temperature  Rating 


Commercial 

Industrial 

Automotive 

Extended 

Military 


0-85  °C 
-40-  100  °C 
-40  -  125  °C 
-40-  125  °C 
-55  -  125  °C 


3.7.4  System  Optimization 

In  parallel  with  the  geometry  optimization  presented  previously,  the  thickness  and  design  of  the 
skin  can  also  be  optimized  for  a  given  objective.  The  above  heating  patterns  and  experimental 
honeycomb  geometry  are  used  for  the  following  demonstrative  case.  Assuming  the  skin  consists 
of  an  unfilled  honeycomb  core  with  filled  honeycomb  on  the  top  and  bottom  surfaces  as  in 
Figure  74,  the  overall  thickness  of  the  skin  and  the  thickness  of  the  SMP  completely  describe  the 
system. 


Figure  74:  SMP/Honeycomb  Skin  Geometry  Showing  Unfilled  Core  and  SMP  Face  Sheets 

A  variational  solution  to  Huber’s  equation  [62],  Equations  (A- 12)  through  (A- 19),  can  then  be 
used  to  solve  for  the  center  point  deflection  of  a  rectangular  panel  with  the  dimensions  and 
boundary  conditions  prescribed  by  the  MAS  program,  Table  1.  The  remaining  MAS 
requirements,  namely  weight,  are  then  calculated  directly  and  implemented  as  constraints. 


D 


X 


d^w 

dx^  dx^dx^ 


+  A, 


(51) 
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(52) 


w  =  ■ 


llq 


2048a^Q 


{4x^  -a^y{Ay^  -a^c^y 


Q  =  25,025(c‘"Z)'  +  D] )+  834,635c^  [c^D,^  +  As )+ 

1 29,740c' (c^Ae  +  Ay )+ 77,363c"  (c"A9  +  Ao)+7,436c^//' 


Ai=7)A 


A2 


Du=D,H 


Du  -D^D^ 


As  =  AA 


Ay  =  A77 


As  =  AA^ 


A9  =  A77^ 


Ao  =  A77^ 


=  a  Wj  +  X  W2  +  j  Wj 


Wi  =1430(c'A  +A)+73,036c"Ai  +  9,477c' (c^Ay  +  As )+ 1183c"//' 
W2  =  1404c"  Ai  +37,180A  +  1144c  "'Ay  +  64,220c' As  +  7436c"//' 


W3  =37,180c^A  +  1404c' Ai  +  64,220c"  Ay  +1144A3  +  7436c'//' 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 


The  optimization  can  then  be  executed  for  each  of  the  heating  patterns  shown  in  Figure  71. 

Figure  75  shows  the  area  density  of  the  optimal  design  for  each  heating  pattern  for  various  total 
thicknesses.  To  generate  Figure  75,  the  total  thickness  is  prescribed  and  the  thickness  of  the 
SMP  layers  is  varied.  The  solid  horizontal  red  line  represents  the  maximum  area  density  as 
prescribed  by  the  MAS  requirements.  If  the  skin  is  required  to  meet  all  MAS  requirements  in 
any  of  the  heating  configurations  except  for  pattern  two,  the  best  possible  design  is  the  thinnest 
total  thickness  where  all  patterns  are  below  the  density  requirement,  occurring  at  a  total  thickness 
of  27.8  mm.  Notice  that  the  curve  for  pattern  two  never  falls  below  the  MAS  density 
requirement  indicating  that  there  is  no  geometry/skin  thickness  meeting  both  the  deflection  and 
weight  requirement  when  the  skin  is  entirely  in  the  soft  state.  Minimum  density  and  skin 
thickness  designs  using  both  PEEK  and  epoxy  are  also  listed  in  Table  20.  As  expected, 
minimum  thickness  designs  result  in  no  empty  honeycomb  while  minimum  density  designs  result 
in  extremely  thick  skins  consisting  of  almost  entirely  empty  honeycomb. 
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Figure  75:  Optimal  Skin  Design  Meeting  MAS  Requirements  for  Each  Heating  Pattern 
with  Respect  to  Area  Density  and  Total  Skin  Thickness 

Table  20:  Optimal  Minimum  Density,  Minimum  Thickness,  and  Best  Overall  Skin  Designs 


Minimum  Density 

Minimum  Thickness 

Best  Overall 

Limiting  Pattern 

3 

3 

3 

SMP  thickness  (m) 

0.00032436 

0.0047197 

0.0014848 

PEEK  Core  thickness  (m) 

0.076249 

4.86E-17 

0.030099 

Max  Deflection  (m) 

0.0025 

0.0025 

0.0025 

Area  Density  (g/m2) 

2561.9942 

12,568.5002 

4624.4238 

Limiting  Pattern 

6 

3 

7 

SMP  thickness  (m) 

0.00065384 

0.0029169 

0.0017489 

Epoxy  Core  thickness  (m) 

0.074139 

5.42E-08 

0.02433 

Max  Deflection  (m) 

0.0025 

0.0025 

0.0025 

Area  Density  (g/m2) 

3179.6866 

6818.85257 

4630.266 

Figure  76  highlights  another  design  using  epoxy  showing  the  center  point  deflection  of  the  panel 
of  each  heating  pattern  with  respect  to  the  MAS  requirement  of  2.5  mm.  All  investigated  heating 
patterns,  except  for  when  the  entire  skin  is  soft,  meet  MAS  requirements. 
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1.7  mm  skin 
33  mm  core 
19,152.1  Pa  Pressure 
0.381  X  0.508  m  Panel 
4633  g/m3  (4638  g/m3  limit) 


- • - ; - 

•  • 

1 ^ ^ ^ ^ _ , 

1  2  3  4  5  6  7  8 

Thermal  Pattern 

Figure  76:  Panel  Center  Point  Deflection  Under  Uniform  Pressure  for  Each  Heating 

Pattern 

3.8  Reconfigurable  Skin  System  Design  Conclusions 

3.8.1  Technological  Advances 

One  of  the  key  benefits  of  the  presented  model  is  its  ability  to  predict  axial  and  shear  stiffness  in 
any  direction.  The  equations  can  easily  be  incorporated  into  an  optimization  scheme  to  take 
advantage  of  the  anisotropies  that  occur  with  honeycomb  geometries  having  different  a  and  / 
member  lengths.  Designs  such  as  those  that  maintain  axial  stiffness  while  minimizing  shear 
strength  in  certain  directions  are  easily  envisioned.  Currently  available  models  calculating  the 
in-plane  stiffness  of  filled  honeycomb  have  been  shown  to  be  inaccurate  at  infill  moduli 
approaching  that  of  epoxy  shape  memory  polymer  above  the  glass  transition  temperature.  The 
presented  experimental  data,  however,  for  filled  honeycomb  above  Tg  correlates  well  with 
moduli  calculated  from  Equations  (38)  and  (42).  This  suggests  that  modeling  the  composite  as  a 
mechanism  where  the  effect  of  the  infill  manifests  as  a  linearly  distributed  load  along 
honeycomb  members  is  appropriate.  While  deformation  due  to  shear  and  torque  and  non-perfect 
geometric  conditions  such  as  rounded  edges  and  non-perfectly  bonded  layers  are  neglected,  the 
presented  equations  are  sufficiently  complex  to  capture  the  physics  of  the  deformation 
experienced  by  the  composite  above  Tg.  As  the  infill  modulus  approaches  that  of  the 
honeycomb,  however,  bending  no  longer  dominates  deformation  and  the  presented  equations 
become  invalid.  In  this  regime,  traditional  rule  of  mixtures  has  been  shown  to  be  sufficiently 
accurate.  To  determine  a  specific  limit,  or  more  likely  to  define  the  accuracy  of  the  model  with 
respect  to  component  relative  modulus,  experimental  data  for  many  more  infill  moduli  is 
required.  With  only  two  data  points,  it  has  only  been  shown  here  that  the  model  is  accurate  at  a 
relative  stiffness  of  3800: 1  and  inaccurate  at  a  relative  stiffness  of  56: 1 .  Also,  the  presented 
linear  model  is  only  valid  for  strains  within  the  linear  elastic  regions  of  the  honeycomb  and  infill 
materials.  From  FEA  analysis  of  a  unit  cell,  the  strain  concentration  ratio  or  maximum  infill 
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strain  relative  to  the  global  strain  was  calculated  to  be  3.4  for  small  deflections,  which  can  be 
used  to  estimate  the  displacement  limitation  of  the  model. 

Also  demonstrated  was  a  novel  heating  scheme  for  shape  memory  polymer.  Aligning  Nickel 
particles  into  chains  using  linear  vibrations  and  a  uniform  magnetic  field  was  shown  to  produce 
relatively  evenly  spaced  conductive  paths  through  the  sample.  Although  untested,  it  is  believed 
that  the  presented  method  could  be  used  for  sample  thicknesses  much  greater  than  the  ones 
shown.  The  method  produced  quick,  even  heat  of  samples  with  minimal  energy  and  minimal 
additional  weight.  Due  to  local  thermal  expansion  of  the  SMP  matrix,  the  nickel  chains  can  also 
be  used  as  an  embedded  temperature  sensor.  Resistance  as  a  function  of  temperature  was 
characterized  and  used  as  a  thermal  switch  indicating  when  a  sufficient  temperature  had  been 
reached  experimentally.  Such  functionality  will  reduce  weight  and  complexity  in  future  systems. 

3.8.2  Future  Work 

Although  many  advances  were  made  during  the  progress  of  this  program,  much  more  research  is 
needed  to  fully  exploit  the  possibilities  of  the  proposed  system.  The  method  used  to  align  the 
nickel  particles  creating  embedded  joule  heating  elements  could  benefit  from  more  detailed 
modeling  and  optimization.  Supporting  technologies  such  as  flexible  electrodes  and  wiring  need 
to  be  investigated  to  properly  integrate  the  system  into  host  vehicle  platforms  and  to  fully  take 
advantage  of  the  strain  capabilities  of  the  skin.  More  research  beyond  what  was  conducted  as 
part  of  this  effort  is  needed  to  increase  the  adhesion  between  the  polymer  and  the  honeycomb, 
which  was  the  main  failure  mode  of  most  samples.  Control  methods  and  algorithms  are  needed 
to  optimally  heat  the  SMP  in  a  manner  consistent  with  desired  motions  of  the  skin.  More  work  is 
needed  to  maintain  heat  where  it  is  needed  preventing  heat  loss  to  the  environment  as  well  as 
spreading  to  areas  of  the  skin  where  it  is  undesired.  Such  thermal  diffiisivity  is  a  driving  factor 
in  the  utility  of  the  skin  and  the  time  available  to  reconfigure  the  structure  under  load.  While 
there  are  issues  that  need  addressed  to  fully  exploit  the  potential  functionality  of  the  skin,  many 
advances  have  been  made  during  the  project  including  the  method  used  to  heat  the  SMP  and  the 
models  and  optimization  schemes  developed  used  to  map  the  performance  of  various  designs 
that  support  further  research  and  implementation  of  the  system. 
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4  PROGRAM  PRODUCTS 


4.1  Publications 

R.V.  Beblo,  JJ.  Joo,  G.W.  Reich, “Joule  Heating  of  Shape  Memory  Polymer  with  Magnetically 
Aligned  Nickel  Particles”,  25’^  International  Conference  on  Adaptive  Structures  and 
Technologies,  Oct  6  -  8  2014,  The  Hague,  The  Netherlands. 

R.V.  Beblo,  J.P.  Puttmarm,N.E.  DeLeon,  J.J.  Joo,  G.W.  Reich,“SMP  Filled  Honeycomb  as  a 
Reconfigurable  Skin:  Model  and  Experimental  Validation”,  19‘^  International  Conference  on 
Composite  Materials,  July  28  -  Aug  2  2013,  Montreal,  Canada. 

R.V.  Beblo,  J.J.  Joo,  G.W.  Reich,  “Aligning  nickel  particles  for  joule  heating  in  epoxy  shape 
memory  polymer  using  a  magnetic  field  and  linear  vibration”  under  consideration  by  the  Journal 
of  Intelligent  Material  Systems  and  Structures. 

J.P.  Puttmann,  R.V.  Beblo,  J.J.  Joo,  B.  Smyers,  G.W.  Reich, “Design  of  a  Morphing  Skin  by 
Optimizing  a  Honeycomb  Structure  with  a  Two  Phase  Material  Infill”,  Smart  Materials, 

Adaptive  Structures  and  Intelligent  Systems  Conference,  Sept  19-21  2012,  Stone  Mountain, 
Georgia,  USA 

R.V.  Beblo,  J.P.  Puttmarm,  J.J.  Joo,  G.W.  Reich,  “Shape  memory  polymer  filled  honeycomb 
model  and  experimental  validation”,  Smart  Materials  and  Structures,  Vol.  24,  No.  2,  January 
2015 

4.2  Patents 

Joo,  James  Jinyong;  Smyers,  Brain  M;  Reich,  Gregory  W;  Beblo,  Richard  Vincent;  Puttmann, 
John  P  “Spatially  targeted  activation  of  a  shape  memory  polymer  based  reconfigurable  skin 
system”  provincial  patent  awarded 
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APPENDIX 


The  forces  due  to  the  infill  when  calculating  the  in-plane  stiffness  of  the  composite  at  angles  a  = 
0,  0+"/2,  and  %  are 
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When  a  is  either  0  or  0+"/2, 11  indicates  the  honeycomb  member  parallel  with  a.  The 
components  contributing  to  deflection  of  the  composite  are  then 
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The  deflection  equation  for  each  of  the  three  directions  becomes 
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Finally,  the  effective  modulus  of  the  composite  for  each  direction  is  then  calculated  using 
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A  similar  process  is  utilized  to  calculate  the  shear  modulus  of  the  eomposite  in  each  direction. 
When  modeling  shear,  the  forces  on  the  honeyeomb  members  due  to  the  infdl  are 
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The  components  contributing  to  the  shear  deflection  of  the  composite  are  then 
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Finally,  the  effective  shear  modulus  of  the  composite  is  calculated  using 
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where  the  components  75  are  solved  for  using  Equations  A-20  through  A-22. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 

Energy-Based  Reconfigurable  Micro-Aerial  Vehicle  (MAV)  Flight  Structures  Design 

A 

Reference  area 

B 

Linearized  measurement  matrix 

Cd 

Coefficient  of  drag 

F, 

Force  along  x  axis 

Gn 

Kalman  gain 

My 

Moment  about  y  axis 

P 

Measurement  error  covariance  matrix 

Qv,w 

Measurement/process  noise  covariance  matrix 

Wb 

Input  weight  matrix 

Wa 

Neuron  connectivity  weight  matrix 

Wb 

Neuron  connectivity  weight  matrix 

Wc 

Output  neurons 

a 

Cylinder  radius 

bn 

Function  denoting  RNN 

n 

Time  step 

V 

velocity  normal  to  hair  sensor 

w 

Vector  of  RNN  connection  weights 

w 

Complex  potential 

a 

Angle  of  incidence  of  the  free  stream 

5 

Drag 

K 

Vortex  strength 

P 

Density  of  air 

Reconfigurahle  Skin  System  Design 

A 

Cross  sectional  area 

Ah 

Area  of  honeycomb  in  a  unit  cell 

Ai 

Area  of  infdl  in  a  unit  cell 

B],  B2,  Di,  D2 

Geometric  parameters 

Cp 

Heat  capacity 

D„Dy 

Flexural  rigidity 

Ecx>  Fey 

Composite  Young’s  modulus 

Eh 

Honeycomb  Young’s  modulus 

Ei 

Infill  Young’s  modulus 

F 

Force 

Fb 

Buoyancy  force 

E, 

Gravitational  force 

Gs 

Honeycomb  shear  modulus 

I 

Second  moment  of  area 

J 

Polar  moment  of  inertia 

M 

Moment 

M 

Magnetic  force 
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N 

Axial  load 

N 

Number  of  cells  in  a  panel 

T 

Temperature 

T 

Torque 

T, 

Glass  transition  temperature 

V 

Shear  force 

V 

Vibrational  force 

Vp 

Particle  volume 

w 

Skin  weight 

a 

Vibrational  amplitude 

a 

Horizontal  wall  length  of  honeycomb 

c 

Depth  of  honeycomb 

d 

Wall  thickness  of  honeycomb 

f(r) 

Initial  temperature  distribution 

k 

Correction  coefficient  for  shear  energy 

k 

Thermal  conductivity 

1 

Slanted  wall  length  of  honeycomb 

m 

Moment 

q 

Applied  pressure 

qin 

Heat  flux 

r 

Particle  radius 

t 

Time 

w 

Panel  center  point  deflection 

a 

Angle  between  unit  cell  orientation  and  x  axis 

5 

Distance  between  two  particles 

5 

deflection 

V 

Poisson’s  ratio 

0 

Interior  half  angle  of  honeycomb 

Pf 

Fluid  density 

Pp 

Particle  density 

T 

Unit  normal  vector 

(0 

Vibrational  frequency 

Abbreviations  and  Acronyms 

AFOSR 

Air  Force  Office  of  Scientific  Research 

AFRL 

Air  Force  Research  Laboratory 

Die 

Digital  Image  Correlation 

EKF 

Extended  Kalman  Filter 

FEA 

Finite  Element  Analysis 

FEMM 

Finite  Element  Method  Magnetics 

IR 

Infrared 

LIDAR 

Laser  Imaging,  Detection,  and  Ranging 

MAS 

Morphing  Aircraft  Structures 

MIMO 

Multi-Input  Multi-Output 
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NI 

National  Instraments 

RNN 

Recurrent  Neural  Network 

SI 

Systeme  Internationale 

SMP 

Shape  Memory  Polymer 

2D 

Two  Dimensional 

UDRI 

University  of  Dayton  Research  Institute 
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